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ATTENUATION  OF  BLAST  WAVES  USING 
FOAM  AND  OTHER  MATERIALS 


INTRODUCTION 


Noise  is  a  problem  everywhere  the  Army  trains  or  tests  with  large  weapons  or 
helicopters.  The  Army  Environmental  Office,  the  U.S.  Army  Construction  Engineering 
Research  Laboratory  (TJSA-CERL),  and  the  U.S.  Army  Environmental  Hygiene  Agency 
(USAEHA)  have  compiled  a  list  of  about  SO  installations  where  training  or  testing  has 
been  altered  or  curtailed  because  of  off -post  noise  problems.  For  these  and  other 
reasons,  the  Army  Environmental  Office  and  the  Secretary  of  the  Army  have  listed  noise 
and  hazardous  waste  as  the  two  major  problems  facing  the  Army  through  the  year  2000 
and  beyond.  Moreover,  these  noise  sources  are  unique  to  the  Army  and  ones  for  which 
technological  solutions  are  for  the  most  part  lacking.  For  U.S.  Army  Europe 
(USAREUR),  noise  is  an  even  bigger  problem  than  in  the  United  States. 

In  general,  USA-CERL  research  is  aimed  at  solving  the  total  noise  problem  for  the 
Army.  This  includes  noise  prediction,  mitigation,  measurement  and  monitoring,  and 
management.  This  report  is  a  compilation  of  articles  on  noise  mitigation,  in  particular 
materials  to  quiet  explosive  noise.  The  articles,  which  have  been  published  in  reviewed 
scientific  Journals  around  the  world,  document  several  years  of  fundamental  research  by 
the  USA-CERL  acoustics  team.  They  are  reprinted  here  in  chronological  order. 
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The  reduction  of  blast  noise  with  aqueous  foam 

Richard  Raspat 

U.S.  Army  Con\trui'tttm  tjiKinrrrmK  Hi'.Si'urch  Lulmnititry,  t\0.  Httx  -ttyty  ( ‘humpmytt,  tlttmu.x  fxtAJif 

S.  K.  Griffiths 

Sandia  National  Laboratories,  P.O.  Box  5800,  Albuquerque,  New  Mexico  87185 

(Received  20  August  1981;  accepted  for  publication  29  August  1983) 

Experiments  were  performed  to  investigate  the  potential  of  water-based  foams  to  reduce  the 
fluffield  noise  levels  produced  by  demolitions  activity.  Measurements  of  the  noise  reductions  in 
.  flat-weighted  sound  etipoaure  level  (FSEL),  C-weighted  sound  exposure  level  (CSBL),  and  peak 
level  were  made  for  a  variety  of  charge  masses,  foam  depths,  and  foam  densities  (230;  I  and  30: 1 
expansion  ratio  foams).  Scaling  laws  were  developed  to  relate  the  foam  depth,  foam  density,  and 
charge  mass  to  noise  reductions.  These  laws  provide  consistent  results  for  reductions  in  the  peak 
level,  FSEL  and  CSEL  up  to  a  dimensionless  foam  depth  of  2.3.  A  two  part  model  for  the 
mechanisms  of  sound  level  reductions  by  foam  Is  suggested. 

FACS  numbers;  43.30.Of,  43.28.Mw 


INTRODUCTION 

To  date,  most  research  into  mitigating  demolition  ef - 
fects  hu  concentrated  on  nearfleld  phenomena,  particularly 
the  damaging  effects  of  the  blasiwave  (the  nearfleld  is  the 
region  within  a  few  hundred  explosive  charge  radii  of  the 
explosive).  Little  hu  been  done  to  meuure  and  evaluate 
blast  elTects  in  the  farfleld,  but  it  is  the  farfleld  elTects  of 
bluts  which  are  becoming  a  serious  environmental  iuue, 
i.e. ,  annoyance  and  damage  complaints  from  individuals  and 
communities  subjected  to  increued  environmental  noise 
levels  can  restrict,  or  eliminate,  blut  producing  activitiu. 

To  address  this  problem  of  the  farfleld  effects  of  demoli¬ 
tion  land  related  activities),  we  have  considered  several 
methods  of  reducing  blast  noise,  including  the  use  of  water- 
based  foams  as  a  mitigating  agent.  This  paper  describes  an 
experimental  investigation  into  the  use  ofboth  low  and  high 
density  aqueous  foams  to  quiet  blast  noise  in  the  farfleld.  The 
data  is  used  to  develop  selling  laws  for  the  foam  so  that  the 
level  of  noise  reduction  can  be  predicted  for  various  amounts 
and  densities  of  foams.  We  also  discuss  past  nearfleld  investi¬ 
gations  when  they  can  be  related  to  our  farfleld  mea.sure- 
ments.  Finally,  a  two  mechanism  model  for  the  reduction  of 
sound  levels  by  foam  is  proposed. 

I.  METHOD 

Two  test  series  were  performed  using  different  expan¬ 
sion  ratio  foams.  The  expansion  ratio  is  the  ratio  of  foam 
volume  to  liquid  content  volume.  To  simplify  the  data  analy¬ 
sis,  simple  lightweight  cubical  or  near  cubical  foam  enclo¬ 
sures  were  used.  For  all  but  the  initial  tests,  these  enclosures 
were  constructed  of  u  wooden  frame  with  polyethelene 
sheeting  for  the  walls.  I'he  charges  were  centered  In  the  cube 
on  crushable  plastic  posts  (Fig.  1).  Spheres  of  Composition 
Four  |C-4)  plastic  explosive  were  used  in  all  tests.  During  all 
tests,  the  charges  were  set  in  pairs:  a  test  charge  under  foam 
and  a  reference  charge  without  foam.  For  cases  where  the 
enclosure  was  slightly  noncubical,  the  geometrically  aver¬ 


aged  foam  depth  was  used  in  the  data  analysis 

</- J(/xuixA)''\  (1) 

where  /.  w,  and  h  are  the  foam  dimensions  in  meters. 

In  all  tests,  the  flat-weighted  sound  exposure  level 
(FSEL),  C-weighted  sound  exposure  level  (CSEL),  and  the 
peak  level  were  measured  at  four  microphone  positions,  two 
each  on  opposite  sides  of  the  explosive.  The  standard  micro¬ 
phone  distances  used  in  most  tests  were  60  and  120  m  (Fig. 
2).  Two  trials  were  performed  for  each  configuration.  The 
levels  were  read  using  a  True  Integrating  Environmental 
Noise  Monitor  and  Sound  Exposure  Level  Meter,  designed 
and  constructed  by  the  U.S.  Army  Construction  Engineer¬ 
ing  Research  Laboratory.  The  signals  were  recorded  for  lat¬ 
er  analysis  on  an  AMPEX  2230  14-track  FM  recorder.  The 
peak  level,  defined  os 

20  lOtip/Pn), 

where /lo  ~  20^  Pa,  is  commonly  used  to  identify  excessive 
noise  levels  around  explosive  facilities  and  to  identify  when 
the  possibility  of  structural  damage  exists.  The  sound  expo- 


FIO.  I.  Experimental  aetup  fur  unuotiflneU  exploalveu  ivai. 
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FIG.  2.  Microphone  poMtionn  for  foum  (esis. 


sure  level  (SEL)  is  deftned  as 

lOlog^ 

where  Ai  deflned  above,  /„  ^  1,0  s  and  the  integral  ia  per* 
formed  over  the  entire  duration  of  the  transient'  The  FSEL  is 
performed  with  no  frequency  weighting,  the  C'Weighted 
SEL  with  a  standard  C-weighting  Alter  in  the  system  before 
the  integration  is  performed. 

The  CSEL  is  important  for  the  environmental  aiseis- 
ment  of  blast  noise  since  it  can  be  used  to  calculate  the  C- 
weighted  average  day/night  level  (CDNL).  The  CDNL  is 
recommended  by  the  Committee  on  Hearing,  Bioacoustics, 
and  Biomechanics  of  the  National  Research  Council'  to  as* 
sess  the  environmental  impacts  of  high-energy  impulsive 
noise  and  la  used  by  the  Department  of  Defense  to  assess 
blut  noise.’ 

II.  SERIES  It  HIGH  EXPANSION  FOAM 

The  flrst  series  of  tests  investigated  the  noise  reducing 
properties  of  high  expansion  foam.  The  foam  was  made  with 
a  National  Foam  System  WP-2S  generator  using  National’s 
l)%  High  Expansion  Foam  solution.  When  water  is  pro¬ 
vided  (by  a  Are  truck)  at  lAQO-ITOO  kPa,  this  generator  pro¬ 
duces  foam  with  a  nominal  250;  1  expansion  ratio.  The  ex¬ 
pansion  ratio  varies  slowly  with  time,  but  samples  taken  30  a 
after  generation  were  usually  within  20%  of  the  nominal 
value.  The  individual  bubbles  were  approximately  I  cm  in 
diameter.  This  type  of  foam  was  stable  and  uuble  for  10-1 S 
min  after  generation  In  low  wind  conditions. 

A.TMt1 

llie  Brst  test  wu  a  feasibility  test.  Two  charge  sizes  and 
two  foam  conAgurations  were  us^.  In  one  case,  the  charges 
were  set  In  a  3.0  m  X  3.0m  x  1.85  m  pit:  the  foam  in  the  pit 
wu  piled  about  0.30  m  above  ground  level.  In  this  partially 
conAned  cau,  the  reductions  In  all  noiu  metrics  were  about 
I4dB  for  the  0.57-kg  charge,  and  about  9  dB  for  the  2.37-kg 
charge. 

The  second  conAguration  wu  an  enclosure 
2.4X2.4X  1.7  m  high,  constructed  u  described  In  Sec.  I. 
This  enclosure  produced  reductions  of  about  10  and  5  dB  in 
all  metrics  (peak  level,  FSEL,  CSEL)  for  the  0.57-  and  2.37- 
kg  chargu,  rupectlvely. 

The  results  of  tut  I  established  that  sIgniAcant  environ¬ 
mental  noiu  reductions  could  be  achieved.  These  reductions 
were  similar  for  all  metrics  measured,  but  did  not  provide 


enough  information  In  allow  foam  thickness  and  charge  sire 
to  he  relaled  to  the  reiiuclion  in  sound  levels, 

B.TMt2 

Test  2  investigated  the  dependence  of  CSEL,  FSEL, 
and  peak  level  on  foam  depth  for  two  charge  masses.  The 
enclosure  dimensions  varied  from  0.30  m  to  1.5  m  in  0.30-m 
steps.  Charges  masses  of  0.57  and  0.061  kg  were  used.  Plots 
of  foam  depth  versus  reduction  for  all  three  metrics  were 
linear  within  the  accuracy  of  the  data.  The  data,  along  with 
those  for  tuts  3  and  4,  for  the  tut  are  displayed  in  Fig.  3 
lA  •«  0.061  kg,  ©  —  0.57  kg)  as  reductions  in  CSEL,  FSEL, 
and  peak  level  versus  the  cube  root  scaled  foam  depth.  The 
results  of  these  six  experiments  show  that  all  three  sound 
levels  are  reduced  linearly  up  to  the  largest  sealed  foam 
depth,  approximately  2.9  m/kg''  *.  At  that  depth,  the  reduc- 


PIO.  .1.  CSEL,  FSEL,  ind  puk  Isvel  rtduuiiuni  veriui  CL-sled  fnaiii  dcpili: 
hl|h  espdiilon.ttllo  foiin  tnli. 
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lions  are  uniformly  about  12  dB.  (In  explosives  work,  it  is 
common  to  use  either  energy  or  mass  scaling  to  display  data. 
In  this  case,  the  scaling  used  is  the  foam  depth  in  meters 
divided  by  the  cube  root  of  the  charge  mass.  This  technique 
will  be  discussed  further  after  the  other  tests  are  described.)’ 

C. TMt3 

Test  3  investigated  the  effect  of  charge  size  fur  a  con¬ 
stant  foam  depth,  The  enclosure  used  In  this  experiment  was 
1 .8  m  on  a  side.  Charge  masses  of 0.061, 0,28, 0.S7, 1. 13,  and 
2.27  kg  were  Ared.  These  data  are  also  displayed  in  Fig.  i 
17  >■  0.28  kg,  a  -  0.37  kg,  -  1.13  kg,  B  =>  2.27  kg).  Be¬ 
cause  the  wind  knocked  the  foam  depth  down  by  0.08  to  0. 1 5 
m  before  the  explosive  was  Ared,  the  geometric  average  of 

l. 7,  1.8,  and  1.8  m  wu  used  In  the  data  analysis.  When  the 
sou.id  level  reductions  are  plotted  versus  the  cube  root 
scaled  foam  depth,  the  plots  increase  linearly  up  to  1.2  m/ 
kg'' at  larger  scaled  foam  depths,  the  reductions  appear  to 
level  off, 

D. T«at4 

When  the  data  from  tests  2  and  3  are  plotted  together,  it 
Is  not  clear  which  results  correctly  describe  the  behavior  of 
the  reduction  at  large  scaled  foam  depths.  The  test  3  data 
show  a  clear  break  in  the  rate  of  reduction  above  10  dB— the 
test  2  data  do  not.  To  investigate  fiirther,  three  experiments 
were  performed. 

ID  A0.06l-kg  charge  was  detonated  in  a  l.6x  l.6x  I.S- 
m  enclosure.  Microphones  were  placed  at  1 3  and  30  m— 
much  closer  to  the  charge  than  in  tests  2  and  3— to  deter¬ 
mine  if  over  land  propagation  differences  could  be  the  cause 
of  the  saturation  and/or  discrepancy  in  results. 

|2)  A  O.S7-kg  charge  was  centered  and  detonated  in  a 
3.7  X  3.7  X  3.7-m  enclosure  with  the  microphones  at  60  and 
120  m. 

(3|  A  2.37-kg  charge  was  centered  and  detonated  in  a 
3. 7x3. 7x3. 7-m  enclosure  with  microphones  at  60  and  120 

m.  Again,  these  results  are  shown  In  Fig.  3, 

E.  Difouaalon  of  toata  2-4 

All  data  from  tests  2-4  are  plotted  In  Fig.  3  versus 
scaled  foam  depth.  The  foam  depths  used  are  the  geometri¬ 
cally  averaged  foam  depths  or  their  equivalents.  The  lines 
shown  are  linear  least  squares  At  to  portions  of  the  data.  The 
Arst  segment  of  the  line  is  Atted  to  the  data  points  from  0.0- 
1.6,  the  second  segment  from  1.2-2. 3.  From  Fig.  3  it  is  ap¬ 
parent  that! 

1 1 1  All  the  data  obey  the  cube  root  scaling  iaw,  with  the 
exception  of  the  single  data  point  from  test  2,  which  lies  well 
above  the  Atted  line.  This  point  is  for  a  small  charge  man, 
0.061  kg,  and  such  charges  are  generally  unreliable.  All  of 
the  other  data  points  are  within  1.3  dS  of  the  lines. 

|2)  The  maximum  possible  reduction  is  limited  to  about 
10  dB.  There  appears  to  be  a  transition  from  a  rapid  increase 
in  reduction  under  1 .2  scaled  m  to  a  much  slower  reduction 
over  1.2  scaled  m.  A  similar  uturation  at  I.S  scaled  m  has 
been  reported  by  Dodly  ti  a/.*  in  their  investigation  of  tl.r 
reduction  of  blast  overpreuures. 

1 7S9  J.  Xooust.  8oo.  Am.,  Vcl.  74,  No.  e,  December  1 063 
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FIO.  4.  Positive  and  loMl  duration  af  a  tran.tleni  waveform  typical  of  an 
expintion  il  d(MC  r«n|c. 


(3)  All  metrics,  CSBL,  FSEL,  and  peak  level,  are  re¬ 
duced  by  roughly  the  same  amount.  The  initial  slopes  are 
3.8,  6.6,  and  6.3  dB  per  scaled  m,  respectively.  The  small 
difference  between  the  peak  level  slope  and  the  FSEL  slope 
Indicates  that  the  foam  causes  little  duration  change  in  the 
wave  form,  since  the  FSEL  is  a  function  of  an  integral  over 
the  duration  of  the  wave. 

To  examine  this  more  closely,  the  positive  phase  dura¬ 
tion  and  total  duration  of  the  wave  were  monitored  directly 
usings  Bruel  and  KJaer  type  7302  Digital  Event  Recorder  to 
display  the  transients  on  a  screen.  The  positive  and  total 
durations  were  measured  from  this  display  {see  Fig.  4).  When 
plotted  versus  scaled  foam  depth,  the  positive  phase  dura¬ 
tion  was  reduced  by  about  3%  by  the  foam;  the  total  dura¬ 
tion  was  reduced  by  about  20%.  Although  there  was  great 
scatter,  the  total  duration  change  tended  to  become  smaller 
as  the  foam  depth  increased. 

To  check  whether  our  farAeld  data  agreed  with  past 
nearAeld  work,  WlnAeld  and  Hill's  data’  were  scaled  out  to 
60  m  using  a  design  chari  of  pressure  versus  distance.  '  This 
technique  Is  at  best  crude,  since  at  close  ronges  energy  is  still 
being  fed  Into  the  shock  wave  by  the  expanding  detonation 
products,  and  the  foam  certainly  must  affect  these  energy 
transfers.  This  calculation  also  neglects  the  reAectlon  of  the 
shock  wave  at  the  foam/air  interface.  Still,  the  results  of  this 
calculation  agree  reasonably  well  with  the  data  from  our 
study.  The  initial  slope  of  the  peak  pressure  reduction  versus 
scaM  foam  depth  line  is  8.6  vs  6. 3  dD  per  scaled  m  measured 
In  the  present  study  .This  difference  may  be  due  to  the  denser 
foam  used  by  WlnAeld  and  Hill,’  or  to  the  placement  of  their 
pressure  transducers  near  the  bottom  of  the  foam  volume, 
where  the  foam  may  be  denser. 

III.  BRnilS  2:  LOW  EXPANSION  FOAM 

A  Mearl  Corporation  OT  10-3  generator  was  used  for 
the  dense  foam  experiments.  The  generator  was  adjusted  to 
produce  a  stiff  30:1  expansion  ratio  foam  at  a  reasonably 
high  Aow  rate;  the  foam  was  made  from  a  3%  solution  of 
National  Foam  Systems  I  high  expansion  foam  solution. 
The  bubble  diameter  in  this  foam  was  on  the  order  of  a  milli¬ 
meter.  I'he  low  expansion  foam  was  quite  stable;  no  drainage 
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or  subtidence  was  noticeable  in  (he  Ant  hour  after  generat¬ 
ing  the  foam. 

The  experimental  set  up  for  the  30: 1  foam  tests  was  the 
same  as  that  described  for  the  high  expansion  foam.  The 
knowledge  gained  from  thr  higli  expansion  foam  results  al¬ 
lowed  a  simpler  experiment  for  the  low  expansion  foam. 

Three  charge  masses  were  used:  0. 1 1 , 0.  S7,  and  2.37  kg. 
Three  cubical  enclosures  were  used  with  the  0. 11  -kg  charge: 
0.3 1, 0.91,  and  1 .3  m  on  a  side.  Five  cubical  enclosures  were 
used  with  the O.S7-kg  charge:  0. 3 1 , 0.6 1 , 0.91 , 1 .2,  and  1 .3  m. 
Two  cubical  enclosure  sizes  were  used  with  the  2.27-kg 
charge:  0.91  and  1.32  m.  The  enclosures  were  oversized  by 
0.2  m  on  length  and  width;  thus,  the  foam  depth  used  in  the 
data  analysis  was  the  geometric  average  of  the  enclosure  di¬ 
mensions  divided  by  two. 

The  reductions  in  noise  level  from  the  various  trials  and 
microphones  were  averaged.  Theae  are  plotted  in  Fig.  3. 


PIO.  3.  CSBt..  PSBL.  ind  peak  lovtl  reduclioni  veriuii  sciled  rosin  depth; 
low  upanilon-rsilo  fotm  tcsii. 


Three  features  in  Fig.  3  are  of  interest. 

(1)  The  data  scale  rather  well;  all  the  points  lie  close  to 
the  best  At  lines  when  reduced  to  scaled  coordinates.  None  of 
jhe  points  lie  further  from  the  Atted  lines  than  1 .0  dB, 

(2/  A”  i-i  the  low  expansion  foam,  the  reduction  in  FSEL 
is  not  linear  over  the  full  range  of  scaled  foam  depth,  but  has 
a  break  point  at  about  0.80  scaled  m.  The  Arst  segment  on  the 
FSEL  curve  is  At  to  the  points  from  0.0-0.9  scaled  m;  the 
second  segment  is  At  to  the  points  from  0.8-1. 3  scaled  m. 
The  CSEL  has  a  similar  break  point  near  0.82  scaled  m.  The 
peak  level  reduction  does  not  display  as  clean  a  break  point; 
however,  the  point  corresponding  to  the  largest  scaled  dis¬ 
tance  is  under  the  curve  At  to  the  rest  of  the  points.  In  each 
case  where  present,  the  break  point  does  occur  at  a  10- 1 1-dB 
reduction — similar  to  the  break  points  in  the  series  I  tests. 

Above  the  break  point,  the  rate  of  reductions  in  the  low 
expansion  ratio  foam  is  greater  than  in  the  high  expansion 
ratio  foam. 

|3)  As  with  the  high  expansion  foam,  the  low  expansion 
foam  reduced  the  FSEL  more  than  it  reduced  the  peak  level. 
The  initial  slope  of  the  peak  level  curve  is  about  11  dB  per 
scaled  m,  while  the  initial  slope  of  the  FSEL  curve  is  about  14 
dB  per  scaled  m.  This  again  indicates  that  the  foam  slightly 
reduces  the  time  duration  of  the  waveforms. 

To  further  investigate  the  characteristics  of  this  reduc¬ 
tion,  positive  and  total  durations  were  measured  in  several  of 
the  tests.  Like  the  high  expansion  foam  test  data,  the  low 
expansion  foam  data  displayed  great  variation  in  duration 
reductions.  The  dense  foam  reduced  the  positive  duration  by 
about  20%:  the  reduction  in  total  duration  was  about  30%, 
with  changes  scattered  down  to  0%  and  up  to  44%.  Even  for 
identical  tests,  the  changes  varied  from  3%-30%.  There  is  a 
small  tendency  for  the  duration  change  to  get  smaller  as  the 
foam  depth  increases,  as  in  the  case  of  high  expansion  foam. 
The  30%  reduction  in  total  duration  corresponds  to  a  1 . 5-dB 
diA'crence  between  peak  and  FSEL  reduction,  if  no  change  in 
shape  occurs. 

IV.  EFFECT  OF  FOAM  DENSITY 

Test  series  I  and  2  considered  only  two  different  foam 
expansion  ratios:  230: 1  (high  expansion  ratio  foam)  and  30: 1 
(low  expansion  ratio  foam).  For  each  foam,  a  cube  root 
scaled  foam  depth  was  used  to  organize  the  test  results  for 
widely  varying  charge  sizes  into  a  single  set  of  curves  for 
each  sound  level  metric.  The  success  in  scaling  the  results  for 
different  chorge  mosses  in  this  way  indicated  that  perliaps 
the  two  seta  of  data  could  be  combined  if  plotted  agaittst  a 
scaled  variable  which  include  the  foam  density.  To  pursue 
this  pcssibility,  the  literature  on  blast  scaling  was  examined. 

An  explosives  scaling  law  which  includes  the  density  of 
the  surrounding  media  is  Lampson's  earth  shock  scaling 
law* 

(/’. -/'o)/Po  =  A(^,/l  Vm),  (2) 

where 

h[  )  is  a  function  only  ofpa(2I  Vw) 
po  is  the  density  of  the  medium  surrounding  the  charge 
R  is  the  distance  from  the  charge  center 
w  is  the  moss  of  the  charge. 


tree  J,  Aoouat.  Soo.  Am,,  Vol.  rx.No.e,  Oeoembenaes  R.  Raspei  and  e  K  Qrinithr  Haduellon  nf  blast  noise  with  foam  l7eo 


9 


Since  the  reduction  data  from  the  teats  of  high  and  low  ex> 
pansion  foamt  scaled  well  as  a  function  of  (/  Vw)*'’,  Lamp- 
son's  scaling  law  suggests  that  plotting  both  seu  of  data  u  a 
function  of  ( p]/w] '  wherep  is  the  foam  density,  may  result 
in  a  uniflcation  of  the  predictions. 

All  of  the  data  from  test  series  I  and  2,  except  for  the 
data  from  teat  1  of  aeries  1,  are  plotted  versus  dimensionleu 
foam  depth  in  Fig.  6.  The  dimensionless  foam  depth  used  in 
this  flgure  is  the  geometrically  averaged  foam  depth  multi¬ 
plied  by  the  cube  root  of  the  iWn  d^ty  in  kilograms  per 
cubic  meter  and  divided  by  the  cube  root  of  the  chvge  mau 
in  kilogiams  of  TNT.  (When  chvge  man  is  used  in  scaling 
laws,  it  is  common  to  expreu  it  in  terms  of  an  equivalent 
mats  ofTNT.  The  C-4us^  in  our  experimenu  is  about  1.34 
times  at  energetic  u  TNT ;  to  agree  with  scaling  conventions, 
our  charge  maatea  were  abutted  by  that  factor.)  The  dimen- 


FIO.  &  CSEL.  FSEL.  and  peak  level  reducdotih  vereuvdimeniionleiii  foam 
depdi. 
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tionless  foam  depth  is 

T  =  (p/Vip)"’,  (3) 

where 

p  is  the  foam  density  in  kg/m’ 

/  it  the  geometrically  averaged  foam  depth 

w  it  the  mats  of  explosive  in  kilograms  of  TNT. 

The  density  of  the  fotim  is  given  by  the  density  of  water  (1000 
kg/m’l  divided  by  the  expansion  ratio. 

From  Fig.  6,  we  tee  that  the  data  scale  well  for  all  me¬ 
trics  up  to  a  dimensionless  foam  depth  oC  2.5.  Little  or  no 
systematic  differences  were  detected  between  the  high  and 
low  expansion  foam  results.  Thus  the  foam  scaling  laws  and 
Pig.  6  can  be  used  to  predict  the  reduction  produced  by  vary¬ 
ing  foam  densities,  foam  depths,  and  charge  masses.  How¬ 
ever,  there  are  not  enough  dau  at  different  foam  densities  to 
establish  that  the  foam  scaling  laws  hold  for  widely  varying 
foam  densities.  For  example,  in  the  extreme  case  of  pure 
water  (expansion  ratio  1:1),  the  foam  scaling  laws  do  not 
hold.  The  reductions  produced  by  water  were  measured  by 
detonating  a  O.S7-kg  charge  of  C-4  in  the  center  of  a  0.39-m 
cube.  The  dimensionless  foam  depth  calculated  for  this  ex¬ 
periment  was  2.22,  which  by  Fig.  6,  would  result  in  reduc¬ 
tions  of  8.0,  S.7,  and  8.2  dB  in  CSEL  and  FSEL,  and  peak 
level,  respectively.  The  actual  average  reductions,  measured 
by  microphones  at  30,  60,  and  120  m  and  averaged,  were 
only  3.8, 3.7,  and  5.7  dB. 

The  foam  scaling  laws  do  not  hold  for  dimensionless 
foam  depths  greater  than  2.5.  Above  2  5,  the  denser  foam 
produced  greater  reductions  than  the  lighter  foam. 

V.  DISCUSSION  OF  RESULTS 

The  results  observed  in  these  tests  can  be  explained  by  a 
two  mechanism  model  of  blut  noise  reduction  by  foams. 
The  noise  reduction  dau  clearly  separates  into  two  regimes, 
above  and  below  a  dimensiuiiieH  foam  depth  of  2.5.  Below 
2.5,  the  data  scales  with  density  and  displays  a  rapid  reduc¬ 
tion  in  noise  level  with  increasing  foam  depth.  Above  2.5  the 
rate  of  reduction  is  smaller  and  does  not  scale  with  density. 
The  break  point  occurs  at  different  scaled  radii  for  the  two 
foams,  however,  it  occurs  at  about  the  same  pressure.  This 
pressure,  calculated  from  the  reduction  daU  and  the  bare 
charge  pressure’  is  on  the  order  of  a  few  hundred  kPa.  From 
independent  tests  of  foam  collapse,  it  has  been  found  that 
these  pressures  are  very  close  to  the  minimum  pressure  nec¬ 
essary  to  fracture  the  foam  cells  and  so  form  a  fine  water 
mist.’  With  this  knowledge  it  is  evident  that  at  leut  two 
distinct  mechanisms  of  noise  reduction  are  operative.  (1) 
Strong  wave  decay  through  the  water  mist  of  the  fractured 
foam;  and  (2)  nonacoustic  decay  of  the  weak  wave  through 
the  intact  foam.  The  first  mechanism  ic  dominant  close  to  the 
explosive,  where  peak  pressures  are  very  high;  the  second 
one  dominates  further  out,  in  the  low  pressure  regime. 

In  general,  the  decay  of  strong  waves  depends  on  the 
irreversible  work  perforinrd  on  the  fluid  between  its  initial 
state  (before  the  wave  has  arrived)  and  its  final  state  (after  the 
fluid  has  returned  to  ambient  pressure)."  ^  This  is  often  called 
the  waste  or  lost  work  and  is  related  to  the  entropy  produced 
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between  (lie  initial  and  final  fluid  states.  For  a  homogeneous 
fluid,  the  source  of  lost  work  is  viscous  dissipation  due  to  the 
very  large  strains  across  the  thin  wave  front.  In  a  two-com¬ 
ponent  material  like  foam,  however,  other  sources  of  irre¬ 
versible  work  are  also  »vailabte;  (he  slip  and  heat  transfer 
between  the  air  and  wuter  as  the  two  components  come  to 
equilibrium  are  both  irreversible  processes.  For  high  peak 
pressures,  the  ^ater  in  the  foam  is  broken  into  extremely  fine 
droplets,  each  with  a  large  specific  surface  area,  so  that  the 
heat  and  momentum  transfer  between  the  two  components 
is  very  rapid.  This  leads  to  substantial  irreversible  work  in 
the  small  time  available  between  (he  arrival  of  the  front  and 
passage  of  the  wave.  This  may  account  for  the  observation 
that  foams  reduce  sound  levels  more  effectively  than  solid/ 
air  systems  of  comparable  density.  That  is,  solids  cannot 
fracture  into  small  particles  that  provide  a  very  large  surface 
area.  In  addition,  the  maximum  possible  strains  are  much 
larger  in  a  foam  than  they  are  in  air  alone;  the  shock  com¬ 
pression  ratio  in  air  Is  limited  to  about  six;  in  dense  foams, 
the  rapid  heat  transfer  from  the  air  to  the  water  permits  a 
maximum  compression  ratio  of  about  30.  Consequently, 
foams  may  exhibit  enhanced  strein-dependent  entropy  pro¬ 
duction. 

Note  that  this  intercomponent  transport  mechanism 
cannot  be  responsible  for  the  observed  reductions  at  large 
foam  depths.  At  the  small  peak  pressures  in  the  foam  at  large 
depths,  the  particle  velocity  and  temperature  rise  are  so 
small  that  heat  and  momentum  transfer  between  the  air  and 
water  cannot  generate  significant  irreversible  work. 

At  large  scaled  foam  depths,  the  peak  pressures  are  very 
small.  In  this  regime,  a  wave  propagating  in  a  homogeneous 
material  would  experience  acoustic  decay  with  (fair  spherical 
geometry)  its  characteristic  inverse  distance  effect  on  the 
sound  levels.  Thus,  if  the  foam  were  homogeneous,  some 
finite  scaled  oepth  would  provide  the  maximum  attainable 
benefit  available  from  a  given  foam.  At  larger  scaled  depths, 
the  sound  levels  would  decay  at  the  same  rate  in  air  and  in 
the  foam  and  so  no  further  reductions  would  be  obtained;  the 
data  in  such  a  case  would  break  to  (brm  a  horizontal  line. 
Since  the  data  show  a  continued  small  increase  in  reductions 
beyond  the  break  point,  we  conclude  that  at  large  scaled 
foam  depths  the  waves  are  experiencing  a  dispersive,  slightly 
nonacoustic  decay  in  (he  additional  foam.  A  weak  wave 
propagating  through  the  foam  is  partially  reflected  at  each 
interface.  This  partial  reflection  results  in  nonacoustic  at¬ 
tenuation  and  dispersion  of  the  wave.  "*  We  will  further  exa¬ 
mine  the  dispersion  when  we  discuss  the  foam  induced 
changes  in  the  pulse  duration. 

We  note  that  the  combined  scaling  of  Fig.  6  supports 
our  two  mechanism  model.  Below  a  dimensionless  foam 
depth  of  2.3  the  data  scale  with  density.  Since  the  foam  is 
shattered  by  the  passing  shock,  the  reduction  depends  only 
on  density  and  not  on  details  of  foam  dimensions  or  cell 
structure.  We  expect  the  scaling  to  hold  in  this  range  pro¬ 
vided  that  the  isothermal  compressibility  of  the  air/water 
mixture  remains  approximately  that  of  air  and  that  direct 
interaction  between  water  particles  (in  the  compressed  state) 
is  not  significant. 

We  likewise  expect  that  the  density  scaling  law  should 


fail  for  large  foam  depths  since  our  proposed  farfield  mecha¬ 
nism  of  low  pressure  weak  decay  depends  not  only  on  foam 
density  but  also  on  foam  structure.  Tlie  low  expansion  foam 
has  a  much  smaller  characteristic  bubble  size  and  so  a  much 
larger  specific  surface  area  than  the  high  expansion  foam. 
Therefore,  it  should  cause  a  more  rapid  attenuation  of  the 
weak  wave  by  reflections  in  the  intact  foam.  This  is  observed 
in  Fig.  3  where  the  attenuation  rate  above  a  dimensionless 
foam  depth  of  2.5  is  much  more  rapid  in  the  low  expansion 
foam. 

The  duration  changes  measured  for  both  foams  also 
.support  our  proposed  two  mechanism  model.  For  the  high 
expansion  foams  the  positive  duration  was  reduced  by  about 
S%  and  the  total  duration  by  about  20%.  The  total  reduc¬ 
tion  becomes  smaller  with  increasing  foam  depths.  For  the 
low  expansion  foam,  the  positive  duration  was  reduced  by 
20%,  and  the  total  duration  reduced  by  30%.  Again,  the 
total  duration  change  becomes  smaller  with  increasing 
depth.  The  large  scatter  in  this  data  prevents  a  detailed  anal¬ 
ysis  of  these  changes,  but  a  consistent  explanation  can  be 
given. 

To  understand  the  reductions  in  wave  duration,  we 
must  address  three  influences:  the  duration  of  the  wave  de¬ 
livered  by  the  explosive,  broadening  of  the  wave  as  it  travels 
through  the  foam  |if  present),  and  broadening  as  it  propa¬ 
gates  through  the  surrounding  air  to  the  microphones.  For 
farfleld  measurements  the  last  of  these — growth  of  the  posi¬ 
tive-phase  duration  as  the  wave  travels  through  the  air — Is 
dominant,  provided  that  the  peak  pressure  at  the  foam/air 
interface  is  not  small."  Most  of  the  positive  duration  is  ac¬ 
crued  in  the  air,  so  that  even  if  the  foam  significantly  altered 
(he  duration  at  the  location  of  the  interface,  this  would  have 
only  a  small  effect  on  the  values  measured  at  a  large  distance 
from  the  charge.  The  small  reductions  noted  in  the  positive- 
phase  duration  are,  therefore,  probably  due  only  to  the  lower 
peak  pressures  (when  foam  is  used)  and  the  consequent  re¬ 
duction  in  rate  of  growth  of  the  pulse  width  in  the  air.  As 
more  foam  is  used,  the  wave  broadening  due  to  the  disper¬ 
sion  inside  the  foam  becomes  significant  enough  to  partially 
offset  thediminishment  due  to  reduced  peak  pressures,  lead¬ 
ing  to  an  overall  smaller  reduction  in  the  positive-phase  du¬ 
ration. 

Since  the  reductions  in  total  duration  are  more  pro¬ 
nounced  than  those  for  the  positive  phase,  it  is  apparent  that 
the  foam  must  have  a  strong  influence  on  the  negative-phase 
duration.  Because  the  negative  duration  is  nearly  indepen¬ 
dent  of  the  length  of  travel  of  the  wave  and  of  the  host  materi¬ 
al,  this  reduction  (unlike  that  of  the  positive  phase)  must 
originate  In  a  reduced  total  duration  delivered  by  the  explo¬ 
sive.  This  is  attributable  to  the  (visible)  absence  of  afterburn 
observed  for  charges  fired  in  the  foam.  Afterburn  occurs 
when  the  detonation  products  are  oxygen  deficient.  These 
products  react  with  the  surrounding  air  to  produce  about 
20%  of  the  total  explosive  energy  of  C-4.  Afterburning  is  a 
relatively  slow  process  and  suppressing  this  reaction  would 
substantially  reduce  the  total  duration.  Although  this  mech¬ 
anism  reduces  the  total  duration  significantly,  the  overall 
energy  reduction  is  small  compared  to  the  reduction  pro¬ 
duced  by  the  other  two  mechanisms  discussed  above.  A  20% 
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direct  reduction  in  explosive  energy  corresponds  to  less  than 
1.0-dB  reduction  in  the  SEL  and  peak  levels. 

VI.  CONCLUSION 

Two  densities  or  aqueous  foam  have  been  tested  for  use 
in  reducing  blast  noise  Trom  a  wide  range  of  explosive 
masses.  The  low  density  250: 1  expansioii>ratio  foam  gives  an 
initial  rate  of  reduction  of  5.8,  6.6,  and  6.3  dB  per  scaled 
meter  of  foam  depth  in  the  CSEL,  FSEL,  and  peak  level 
measurements,  respectively.  In  each  sound  metric,  a  break 
occurs  in  the  reductions  at  about  10  dB;  this  takes  place  at  a 
scaled  foam  depth  of  1.5  m/kg'^^  and  additional  foam  be¬ 
yond  this  point  gives  a  much  smaller  rate  of  noise  reduction. 

The  high  density  30: 1  expansion-ratio  foam  produces 
initial  rates  of  reduction  of  14.0,  IS.I,and  1 0.9  dB  per  scaled 
meter  of  foam  depth  in  the  CSEL,  FSEL,  and  peak  sound 
levels.  As  with  the  low  density  foam  a  break  occurs  in  the 
rale  of  reduction  of  the  FSEL  and  CSEL  at  a  10-dB  total 
reduction.  In  the  high  density  case,  however,  the  break  takes 
place  at  a  smaller  scaled  depth  of  about  0.9  m/kg' No 
clear  break  in  the  rate  of  reduction  of  the  peak  level  was 
found  for  the  high  density  foam. 

When  all  the  test  data  are  scaled  by  the  foam  density  as 
well  as  the  explosive  mass,  the  reduction  results  for  both  the 
high  and  low  density  foams  fall  on  a  single  curve  for  each 
sound  metric— up  to  a  dimensionless  foam  depth  of  2.5.  Be¬ 
low  this  dimensionless  foam  depth,  our  proposed  mecha¬ 
nism  of  noise  reduction  is  the  strong-wave  decay  due  to  irre¬ 
versible  intercomponent  heat  and  momentum  transfer 
between  the  air  and  water  in  the  foam.  This  mechanism  de¬ 
pends  only  on  the  density  and  not  details  of  the  foam  struc¬ 
ture. 

Above  a  dimensionless  foam  depth  of  2.5  the  high  den¬ 


sity  foam  gives  a  larger  rate  of  reduction  than  the  low  density 
foam.  In  this  regime,  the  apparent  mechanism  of  mitigation 
is  the  reflection  of  the  waves  within  the  intact  foam.  The  high 
density  foam  gives  a  larger  rate  of  reduction  here  only  L)e- 
cause  of  its  smaller  cells  and  larger  internal  surface  area.  The 
rates  of  reduction  in  the  sound  metrics  is  much  smaller  in 
both  foams  above  a  foam  depth  of  2.5,  indicating  that  the 
strong-wave  mechanism  is  significantly  more  important 
than  the  tionacoustic  decay  mechanism  for  noise  reduction. 
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Summary 

Noise  produced  by  biast  waves  can  be  a  probiem,  especlaliy  when  an  explosion  occurs 
near  populated  areas.  As  one  means  of  reducing  the  blast  noise,  the  explosive  is  detonated 
in  a  pit,  a  space  closed  at  the  bottom  and  sides,  open  at  the  top.  A  two-dimentional  finite 
difference  model  was  used  to  simulate  such  an  explosion  in  a  pit  and  to  determine  to 
what  extent  the  blast  wave  was  attenuated.  The  code  used,  CSQ,  developed  by  scientists 
at  Sandia  National  Laboratories,  was  tailored  for  our  studies,  '^e  key  results  were:  (a) 
the  presence  of  a  pit  in  all  cases  caused  the  blast  wave  to  be  attenuated;  (b)  for  a  cy> 
lindrical  pit,  a  pit  of  a  radius  which  effected  maximum  blast  wave  attenuation  was 
found;  (c)  a  useful  parameter,  d£out/dfi  energy  loss  rate  from  the  pit,  was  shown  to 
be  a  good  indicator  of  relative  pit  effectiveness. 


Introduction 

It  is  possible  to  reduce  the  noise  and  the  possible  blast  damage  resulting 
from  an  explosive  reaction  by  placing  the  explosive  material  in  a  partially 
open,  thick-walled  container  or,  in  other  words,  in  a  pit.  Examples  of  appli¬ 
cations  would  include  finding  ways  to  reduce  the  blast  noise  near  blasting 
sites  used  by  geologists  to  find  underground  oil  and  gas,  and  determining 
appropriate  ways  to  dispose  of  explosive  armaments  near  populated  areas. 

Questions  that  need  to  be  answered  to  determine  ways  to  utilize  pits  sur¬ 
rounding  an  explosive  to  mitigate  blast  effects  include  the  following: 

(1)  What  is  the  optimum  horizontal  distance  from  the  charge  to  the  pit  walls 
in  order  to  minimize  the  leading  edge  pressure  at  a  fixed  observer  loca¬ 
tion? 

(2)  What  is  the  minimum  vertical  distance  (floor  to  open  top)  to  effectively 
reduce  the  blast  noise  at  a  fixed  observer  distance  from  the  pit? 

(3)  What  parameters  (other  than  the  peak  pressure  pulse  measured  at  an 
observer  location)  can  be  calculated  to  design  the  most  effective  pit? 

(4)  What  is  the  effect  of  energy -absorbing  pit  wstlls? 
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In  this  study  w«  modelled  the  dynamic  processes  in  and  surrounding  the 
pit  by  providing  the  solution  to  the  time-dependent  equations  which  con¬ 
serve  mass,  momentum,  and  energy  of  the  gases  produced  by  the  explosive 
as  well  as  the  air  inside  the  pit.  Clearly,  a  minimum  of  two  independent 
space  dimensions  is  required.  The  blast  produces  shocks  and  strong  compres¬ 
sion  waves  which  propagate  at  supersonic  speeds  relative  to  the  undisturbed 
air.  The  waves  can  move  several  kilometers  per  second.  A  km/s  is  also  one 
mm/jjis,  and  therefore  it  follows  that  the  flow  variables  pressure,  density, 
and  energy  need  to  be  calculated  at  microsecond  (10~‘  s)  time  levels  since 
changes  in  the  flow  variables  occur  over  distances  in  the  scale  of  millimeters 
or  less. 

The  code  CSQ,  developed  by  scientists  at  Sandia  National  Laboratories, 
documented  in  a  Sandia  Technical  Report  [1],  was  used  to  model  the  prob¬ 
lem. 


Governing  equations 

The  two-dimensional,  time-dependent,  Lagrangian  form  of  the  governing 
equations  of  conservation  of  mass,  momentum,  and  energy  in  cylindrical 
coordinates  used  are 


Conservation  of  mass 


Itu-  ri  +  iiiil 
at  ""^Lr  dr  dzJ 

Conservation  of  momenta 
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Conservation  of  energy 
dE  ~  a  /  1  \ 

,2,3) 


Here  p  represents  density:  t,  time;  r,  the  radial  coordinate;  Vr,  the  particle 
velocity  in  the  x  direction;  z,  the  height  coordinate;  Vg,  the  particle  velocity 
in  the  z  direction;  d/dr  and  d/dz,  the  partial  derivatives  with  respect  to  the 
Lagrangian  coordinates,  r  and  *,  respectively:  d/dt,  the  partial  derivative  with 
respect  to  time;  o,  the  stress  tensor;  Q,  the  artificial -viscosity;  P,  the  material 
pressure;  and  E,  the  specific  material  energy.  These  equations  must  be  sup¬ 
plemented  with  an  equation  of  state  to  have  a  determinate  system, 
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Applicution  of  CSQ  to  the  ezploiion*ln-pit  pfoblem 

Am  fteted  ettrlier,  the  program  CSQ  hea  been  used  to  model  the  detonation 
of  a  high  explosive  aurrounded  by  air  in  a  cylindrical  pit  open  to  the  atmo¬ 
sphere.  Figure  1,  the  half-iection  view  of  a  right  circular  cylinder,  shows  the 
features  and  gtidding  of  a  typical  configuration.  By  placing  the  explosives  at 
the  center  of  the  pit,  it  is  possible  to  reduce  the  problem  from  three  dimen¬ 
sions  to  two,  with  symmetry  in  the  angular  {9 )  direction. 


PtB-  1.  Dstcrlptlons  of  th«  grid  uNd  by  CSQ  to  inodtl  otMt  of  varying  oharga  maaa,  am* 
biant  tamparatura,  and  pit  radlua. 

Two  types  of  boundary  conditions  were  employed,  reflecting  and  trans¬ 
mitting.  In  all  cases  considered  in  this  study,  pit  walls  were  treated  u  per¬ 
fectly  reflecting  boundaries.  Velocities  at  reflecting  boundaries  were  defined 
to  be  zero.  At  transmitting  boundaries,  mass,  monientum,  and  energy  fluxes 
were  permitted.  CSQ  used  appropriate  extrapolation  schemes  to  determine 
the  magnitudes  of  the  fluxes. 

CSQ  is  a  two-dimensional  (planar  or  cyllndrlcsil)  hydrodynamic,  Lagrang- 
lan  code  which  Is  capable  of  solving  many  problems  given  an  equation 
of  state  for  each  material  and  proper  initial  and  boundary  conditions.  Given 
a  configuration,  CSQ  will  use  finite  difference  techniques  to  advance  the 
system  in  time.  Details  of  the  differencing  scheme  employed  by  CSQ  are 
found  on  pp.  Td—S?  of  Ref.  [1],  Acceleration,  velocity,  position,  density, 
and  cell  volume  of  each  Lagrangian  cell  are  datermined  using  explicit  finite 
difference  relations.  The  energy  equation  and  equation  of  state  are  then 
solved  simultaneously  using  what  'a  essentially  i  Newton— Raphson  tech¬ 
nique.  At  this  point  all  thermodynamic  quantities  are  known. 
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In  order  to  run  CSQ,  the  user  mutt  define  the  problem.  This  is  accom¬ 
plished  by  executing  the  program  CSQGEN,  in  which  the  geometry,  grid- 
ding,  flow  variables,  thermodynamic  properties  and  equation  of  state  infor¬ 
mation  for  each  material  are  defined.  The  finite  difference  grid  is  fixed.  It  is 
allowed  to  move  during  Laigrangian  calculations  but  it  then  rezoned  back  to 
its  original  position.  This  rezonc  gives  the  code  an  Eulerlan  nature. 

CSQQEN  alto  requires  equation  of  state  information  for  each  material, 
tn  the  case  studied,  this  meant  providing  tabular  data  for  air  and  appropriate 
constants  for  the  JWL  [2]  form  of  the  equation  of  state  for  the  explosive, 
TNT.  Two  thermodynamic  properties  of  each  material  mutt  be  defined  as  an 
Initial  condition  in  CSQOEN  along  with  the  velocity  of  each  material. 

The  tabular  data  for  air  were  provided  by  Sandia  National  Laboratories. 
The  data  describe  the  nonideal  nature  of  air  at  high  temperatures.  Data  were 
available  for  temperatures  ranging  from  232  K  to  2.82  X  10‘  Kand  for  densi¬ 
ties  ranging  from  IX  10' *  g/cm’  to  3  X  10~^  g/cm\  This  yielded  data  for  a 
range  of  pressures  from  1.0  dyne/cm*  to  1  X  10'*  dyne/cm*.  For  all  prob- 
;  'ms  studied,  the  properties  of  the  air  remained  within  these  limits. 

Among  the  many  quantities  that  were  varied  in  the  system  studied  were 
pit  geometry,  charge  mass,  and  ambient  temperature  of  the  surrounding  gas 
[3].  However,  each  problem  studied  retained  certain  basic  features,  namely, 

(1)  A  spherical  charge  of  a  high  explosive  sharing  a  common  centerline  with 
a  right  circular  cylinder,  lying  on  the  base  of  the  cylinder  was  detonated, 

(2)  Interactions  occurred  when  the  spherical  wave  generated  by  the  ex¬ 
plosion  struck  a  reflecting  boundary,  and 

(3)  Fluid  motion  occurred  outside  the  pit. 

Detonation 

The  detonation  of  a  high  explosive  is  characterized  by  the  progression  of 
a  detonation  wave  which  moves  at  a  constant  velocity,  D,  and  which  leaves 
the  reacted  material  at  a  pressure,  Poj»  the  Chapman— Jouguet  pressure.  The 
variables  D  and  Pqj  are  specific  to  the  explosive  and  are,  in  general,  calcula¬ 
ble  quantities.  Details  on  this  are  well-documented  [4]  and  will  not  be 
repeated  here.  The  propagation  of  the  detonation  wave  through  the  high 
explosive  was  modelled.  The  behavior  of  a  spherical  blast  wave  which  is 
surrounded  only  by  air  is  well-known.  The  problem  is  one-dimensional  in 
the  radial  direction  in  spherical  coordinates.  Relations  which  predict  pres¬ 
sure  as  a  function  of  radius  are  given,  for  example,  by  Baker  [5]. 

CSQ  was  used  to  simulate  the  unconfined  explosion  of  a  spherical  charge 
of  TNT,  r  1  cm,  in  air  at  standard  atmospheric  conditions.  Since  the  code 
is  a  finite  difference  code  in  cylindrical  coordinates,  a  true  sphere  cannot  be 
defined.  The  spherical  charge  is  simulated  by  a  set  of  cylindrical  finite  dif¬ 
ference  cells.  Cells  of  mixed  composition  are  allowed,  so  it  was  possible  to 
have  the  mass  of  explosive  which  corresponded  to  the  specified  charge 
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radius  and  density.  Cells  of  mixed  composition  were  treated  as  a  homo¬ 
geneous  mixture;  no  internal  boundaries  existed. 

Results  of  this  calculation  are  shown  in  Fig.  2,  along  with  three  results 
obtained  by  other  means.  The  peak  pressure  at  a  given  distance  from  the 
charge  center  is  plotted  as  a  function  of  the  number  of  charge  radii  away 
from  the  charge  center.  CSQ*s  predictions  were  arbitrarily  made  (dong  the 
vertical  centerline.  Two  of  the  other  predictions  shown,  found  in  Baker’s 
book  [5],  are  based  on  experiments.  The  other  predictions  shown,  made 
by  Griffiths  [6],  is  based  on  calculations  that  assumed  a  radially  one-di¬ 
mensional  pressure  wave.  It  is  seen  that  none  of  the  four  methods  agree 
perfectly,  but  that  CSQ  predicts  pressures  in  the  same  range  as  those  pre¬ 
dicted  by  the  one-dimensional  methods. 


Pig.  2.  Predictions  (or  the  preMure  st  the  leading  edge  of  the  rsdisliy  one-dimensional 
spherical  blast  wave  by  various  methods. 


InteracUon  of  waves  with  boundaries 

As  stated  earlier,  complex  Interactions  between  spherical  shock  waves  and 
reflecting  surfaces,  both  planar  and  cylindrical,  occurred  in  the  problem 
modelled.  In  Ref.  [3]  the  CSQ  code  was  used  to  model  a  simple  case  of  a 
shock  wave  interacting  with  a  reflecting  botmdary.  The  problem  modelled 
was  the  classical  one-dimensional,  linear  shock  tube,  divided  into  two  re¬ 
gions,  one  of  high  pressure  air  and  the  other  of  low  pressure  air.  The  system 
was  assumed  to  be  in  thermal  equilibrium.  The  results  obtained  from  the 
solution  of  the  analytic  equations  given  by  Shapiro  [7]  and  the  results 
from  the  application  of  CSQ  were  almost  identical  for  the  pressure  behind 
the  initial  shock  and  the  pressure  behind  the  shock  reflected  from  the 
closed  end.  Since  the  linear,  one-dimensional  reflections  were  accurately 
modelled,  it  was  assumed  that  CSQ  could  accurately  predict  more  com¬ 
plicated  reflections  such  as  those  that  result  when  explosions  occtu  in  pits. 
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Complicated  two'dimensional  wave  interactions  were  modelled.  The  wave 
motion  is  easily  described  in  the  early  stages  by  a  simple  spherical  blast 
wave.  Upon  striking  the  cylindrical  wall,  a  shock  wave  is  reflected  towards 
the  pit  center.  This  shock  interacts  with  the  initial  blast  wave.  Soon  after 
this,  all  combinations  of  interactions  (shock-shock,  shock-expansion, 
shock*-wall,  etc.)  are  possible.  An  accurate  verbal  description  becomes 
nearly  impossible. 

Two  limiting  cases  can  be  easily  described.  With  JR  as  the  radius  of  the 
pit  and  JFf  as  its  height,  the  parameter  A/H  is  a  characteristic  of  the  pit.  As 
RjH  approaches  infinity,  an  explosion  in  the  pit  behaves  essentially  as  if 
there  were  no  pit;  it  is  really  an  unconfined  explosion  except  for  the  inter¬ 
action  with  the  horizontal  surface,  which  only  serves  to  double  the  effective 
charge  energy.  At  the  other  end,  as  RIH  approaches  zero,  the  blast  wave 
attains  a  linear  nature.  Reflections  occur,  and  the  resulting  wave  looks  like 
a  planar  wave  propagating  through  a  tube. 

As  the  blast  wave  leaves  the  pit,  it  begins  to  behave  spherically.  Behind  the 
initial  front,  waves  ore  still  interacting,  and  they  too  leave  the  pit.  In  some 
cases  studied  it  was  observed  that  these  waves  catch  up  to  and  strengthen  the 
initial  blast  wave.  A  feature  noticed  in  all  cases  was  a  vortex  which  developed 
on  the  surface  outside  the  pit  near  the  pit  wall. 

A  final  feature  of  problems  studied  is  the  relatively  long  computation  time 
required  to  model  the  problem,  In  most  cases  the  maximum  number  of  cells 
was  used  in  order  to  achieve  the  most  accurate  results.  As  more  cells  are 
used,  computation  time  to  model  equivalent  systems  increases.  To  model  the 
detonation  of  a  charge  of  TNT  required  about  1000  CP  seconds  on  the  Uni¬ 
versity  of  Illinois’s  CYBER  175  computer.  Once  this  foundation  was  built, 
tests  varying  pit  geometry  could  be  made. 

Results 

Numerical  tests  were  conducted  which  predicted  the  history  of  the  cumu¬ 
lative  amount  of  energy  which  had  left  the  pit.  Other  numerical  tests  exam¬ 
ined  the  flow-field  both  in  and  outside  of  the  pit.  The  parameters  which 
were  varied  were  (a)  pit  height,  (b)  pit  radius,  (c)  charge  mass,  and  (d) 
ambient  temperature. 

For  this  study,  we  define  an  effective  pit  on  a  relative  scale.  For  the 
general  problem  of  the  detonation  of  high  explosives  within  a  partially  en¬ 
closed  volume,  one  can  state  that  the  more  the  pressure  wave  that  results 
from  the  detonation  is  attenuated,  the  more  effective  is  the  partially  en¬ 
closed  volume.  As  stated  before,  the  partially  enclosed  volumes  studied 
were  cylindrical  volumes  which  were  closed  at  the  base,  closed  at  the  sides, 
and  open  at  the  top;  or,  more  compactly,  pits.  The  tests  were  limited  to 
models  of  explosions  of  spherical  charges  of  TNT  lying  on  the  bottom  of  the 
cylindrical  volume  at  its  centerline. 
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Flow-fitid  inside  pit:  Eout  versus  time 
As  one  means  of  determining  pit  effectiveness,  the  time  derivative  dffout/ 
df  is  calculated.  £out  i*  the  cumulative  amount  of  energy,  both  internal  and 
kinetic,  which  has  left  the  pit.  This  number  results  from  the  utilization  of 
the  transmitting  boundary  condition  defined  to  exist  at  the  top  of  the 
mesh.  The  transmitting  boundary  condition  allows  fluxes  in  and  out  of  the 
mesh  by  defining  flow  variables  outside  of  the  mesh  to  be  extrapolated 
values  of  the  flow  variables  in  the  mesh. 

We  believe  that  as  pit  effectiveness  increases,  energy  release  decreases, 
causing  the  downstream  pressure  pulse  to  be  weaker.  A  less  effective  pit 
would  allow  a  rapid  release,  causing  the  pressure  pulse  to  more  closely 
resemble  the  pressure  pulse  that  would  exist  had  no  pit  been  present.  The 
measure  of  pit  effectiveness  dfout/df  1*  ^  qualitative  measure.  It  is  useful  for 
cases  in  which  a  single  parameter  is  varied. 

Eout  versus  time  for  explosions  in  pits 
Three  pit  heights  were  modelled  in  these  tests.  Pit  radius  was  main¬ 
tained  at  100  cm.  A  charge  of  TNT,  r  ■  1  cm  (mass  ■  6.83  g),  was  used  in 
all  cases,  and  ambient  conditions  of  the  sunounding  air  were  taken  as 
standard  atmospheric.  The  three  pit  heights  modelled  were  0.4, 1.2,  and  2.0 
m.  To  eliminate  any  possible  grid  biasing  between  the  three  cases,  a  constant 
cell  size  was  used,  namely  6.714  cm  by  2.857  cm.  This  resulted  in  variable 
grid  dimension  for  each  case,  36  X  36  for  the  2.0  m  height,  21  X  36  for  1.2 
m,  and  7  X  35  for  0.4  m. 


Fis.  8.  esq  prtdiotion  o(  vtriui  tiin<  for  datonatlona  in  oyllndrloal  pits  of  variabit 
height. 
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The  predictions  for  £out  are  shown  in  Fig.  3.  Here,  is  plotted  as 
a  function  of  (t  —  to),  where  to  is  the  time  when  energy  was  fmst  predicted 
to  escape  the  pit. 

As  the  results  clearly  indicate,  dEout/dt  increases  as  pit  height  decreases. 
For  any  given  time,  more  energy  will  have  left  the  shorter  pit  than  the  taller 
pit.  From  this  we  conclude  that  if  a  blast  wave  is  partiaily  enclosed  by  a 
cylindrical  pit  with  reflecting  walls,  as  the  wall  height  increases,  the  blast 
wave  appears  to  have  been  produced  by  a  charge  of  smaller  mass.  In  sum¬ 
mary,  as  pit  height  incieases,  dEout/df  decreases,  and  pit  effectiveness  in¬ 
creases. 

Flow-fitld  outsidt  pit:  Pu  versus  r 

In  order  to  obtain  qualitative  results,  namely,  the  flow-field  variables 
outside  the  pit,  it  was  decided  to  expand  the  domain  of  the  problem  to  in¬ 
clude  the  volume  outside  the  pit.  In  all  further  tests,  a  constant  grid  geom¬ 
etry,  shown  earlier  in  Fig.  1,  was  maintaiiied.  The  number  of  cells  used, 
4000,  corresponding  to  a  50  X  80  grid,  is  close  to  the  maximum  number 
available  in  our  CYBER  176  computer. 

Variable  pit  radius 

A  series  of  tests  was  performed  that  modelled  the  detonation  of  a  sphere 
of  TNT,  r  ■  1  cm,  surrounded  by  air  at  standard  atmospheric  conditions.  The 
charge  was  located  at  the  bottom  center  of  a  cylindrical  pit  of  constant 
height,  100  cm,  and  variable  radius.  Sbc  cases  were  modelled,  with  the  pit 
radii  in  the  tests  being  20,  26,  36,  60, 76,  and  160  cm.  The  geometry  of  the 
configuration  is  teen  in  Fig.  4. 

Figure  6  shows  tome  pressure— distance  results  of  these  tests.  Here  the 
predicted  pressure  of  the  leading  edge  of  the  wave  at  points  on  the  hori¬ 
zontal  surface  outside  the  pit  is  plotted  as  a  function  of  distance  from  the 
bottom  center  of  the  pit  for  pits  of  various  radii.  Also  shown  are  results 
for  r  •>  infinity,  which  corresponds  to  the  case  where  no  pit  is  present; 
only  the  horizontal  surface  exists.  Far  from  the  charge,  this  case  is  equiva¬ 
lent  to  the  unconflned  explosion  of  a  charge  of  twice  the  original  charge’s 
ma«s. 

By  the  previous  definition  of  an  effective  pit,  we  say  that  at  a  given 
distance  from  the  charge,  the  lower  the  pressure,  the  more  effective  the  pit. 

Some  key  features  illustrated  on  this  figure  are: 

(1)  the  decay  of  the  leading  edge  pressure  with  increasing  distance  from  the 
charge  center  is  predicted  in  all  cases  studied, 

(2)  at  all  points  on  the  horizontal  surface  outside  the  pit,  the  existence  of 
vertical  barrier  walls  is  predicted  to  make  the  pit  more  effective  than  a 
pit  with  no  barrier  walls,  and 

(3)  the  existence  of  a  minimum  leading  edge  pressure  as  a  function  of  pit 
radius  is  predicted.  That  is,  it  is  predicted  that  at  a  given  distance  R 
from  the  charge  center,  the  derivative  dPb/drpit  is  zero  when  evaluated  at 
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a  certain  barrier  radius,  and  this  critical  point  is  a  point  of  minimum 
pressure. 

The  existence  of  the  minimum  is  seen  more  clearly  in  Fig.  6  in  which  the 
pressure  at  the  leading  edge  of  the  wave  at  two  points  along  the  horizontal 
surface  (outside  the  pit)  is  plotted  as  a  function  of  pit  radius.  At  this  dis* 
tance,  the  most  effective  pit  has  a  radius  of  approximately  50  cm. 


Pig.  6.  Prodictsd  praMurs  «t  leading  adga  at  i  point  on  tha  lurfaoa  outilda  the  pit  variui 
barrier  rtdiui. 


Point  3  deserves  some  amplification  and  speculation.  First,  the  existence 
of  this  minimum  was  not  expected.  What  was  expected  was  that  pit  effective" 
ness  would  increase  as  pit  radius  decreased.  A  pit  of  infinite  radius  is  really 
no  pit  at  all;  it  is  an  unconflned  region.  As  the  radius  is  brought  in  from  in¬ 
finity  to  zero,  it  was  expected  that  the  pit  would  have  more  and  more  of  a 
damping  effect  on  the  blast  wave.  The  increased  damping  effect  was  pre¬ 
dicted  by  CSQ,  but  only  to  a  point  at  which  other  factors  must  have  in¬ 
fluenced  the  flow’s  behavior. 

One  possible  explanation  is  that  the  original  expectations  are  indeed  true 
in  the  far-fleld.  Possibly  the  domain  of  the  problem  studied  is  not  large 
enough  to  observe  this,  and  therefore  it  may  be  true  that  only  a  near-field 
phenomenon  has  been  observed. 

Another  possible  explanation  of  the  minimum  is  related  to  secondary 
waves  trailing  the  leading  edge  wave.  In  all  cases,  CSQ  describes  secondary 
waves  reflected  from  the  various  surfaces  in  the  domain.  It  is  noted  in  some 
cases  that  these  secondary  waves  leave  tne  fluid  at  a  higher  pressure  than  the 
leading  edge  wave  did  at  the  same  point.  As  pit  radius  shrinks,  it  is  predicted 
that  the  secondary  waves  overtake  and  stre>igthen  the  leading  edge  wave. 
This  presents  the  possibility  that  in  the  far-field,  all  secondary  waves  will 
eventually  strengthen  the  leading  edge  wave,  possibly  rendering  any  attenua- 
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tion  by  the  pit  to  be  nonexistent.  This  is  plausible  for  the  system  studied, 
since  the  reflecting  walls  do  not  remove  any  energy  flrom  the  blast  wave. 

Regardless  of  why  the  minimum  is  observed,  a  consequence  of  its  exis* 
tence  is  important.  Because  of  the  minimum,  if  one  were  designing  a  pit  to 
attenuate  a  blast  wave,  an  optimum  pit  radius  could  be  determined. 

Variable  eharg*  man 

Simulations  testing  the  effect  of  charge  mass  were  carried  out.  In  these 
tests,  the  detonation  of  a  spherical  charge  of  TNT,  r  ■>  8  cm  (mass  ■  3497  g), 
rather  than  r  1  cm,  was  modelled.  Again,  the  charge  was  surrounded  by 
standard  atmospheric  air  and  a  pit  of  constant  heighti  h  -  100  cm,  and 
variable  radius,  r  ■  20, 26,  35,  50, 100, 160  cm. 

Figure  7,  similar  to  Fig.  5,  shows  the  predicted  pressure  at  the  leading 
edge  of  the  wave  at  points  along  the  horizontal  surface  outside  ^e  pit  as  a 
function  of  distance  from  the  bottom  center  of  the  pit  for  pits  of  various 
radii,  including  infinite  radius.  The  same  features  noticed  in  the  earlier 
series  of  tests  (with  a  smaller  explosive  charge),  decay  of  the  leading  edge 
pressure  with  distance  from  charge  center,  the  ability  of  the  pit  to  attenuate 
the  blast  wave,  and  the  existence  of  a  barrier  radius  which  yields  a  minimum 
leading  edge  pressure,  are  all  predicted  hen.  However,  at  these  higher  blast 
preuuns  one  notices  that  pndicted  behavior  is  leu  explainable.  The  pres¬ 
sure-distance  results  exhibit  different  featuns  than  the  results  obtained  for 
the  explosion  of  the  smaller,  r  ■  1  cm  charge.  For  example,  a  point  of 
maximum  pnuure  at  a  given  barrier  radius  is  predicted  for  rpu  <  36  cm.  The 
only  explanation  offered  for  this  is  that  stronger  shock  waves  tend  to  make 
systems  leu  linear,  and  thus,  generally  much  more  difficult  to  describe. 


Fig.  7.  Predicted  prewure  on  surface  outside  of  pit  of  the  leading  edge  versus  distance 
from  the  bottom  center  of  the  pit  for  pits  of  variable  radius. 
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Variable  ambient  temperature 

Another  test  conducted  varied  the  ambient  temperature  of  the  surround¬ 
ing  air.  Three  cases  were  studied,  one  with  the  temperature  of  a  cold  day, 
246  K  ('-18.4*’F),  one  at  standard  atmospheric  conditions,  298  K  (77.0‘’F), 
and  one  with  the  temperature  of  a  hot  day,  320  K  (116.6°?).  In  all  cases, 
the  pressure  was  maintained  at  standard  atmospheric  level  and  the  density 
adjusted  to  satisfy  the  equation  of  state  for  air.  The  detonation  of  a  charge 
of  TNT,  r  "  1  cm,  lying  at  the  bottom  center  of  a  pit  radius  25  cm,  height 
100  cm  Is  modelled.  In  all  three  cases,  it  was  seen  that  ambient  temperature 
had  no  effect  on  the  pressure  at  the  leading  edge  of  the  wave  in  the  domain 
studied. 

Conclusions  and  Recommendations 

The  main  conclusions  that  can  be  drawn  from  the  results  of  this  study  are; 

(1)  By  using  the  two<dlmensional,  axisymmetric  finite  difference  form  of  the 
governing  equations  of  conservation  of  mass,  momentum,  and  energy, 
along  with  equations  of  state  for  the  materials  studied,  it  is  possible  to 
model  explosions  and  describe  interactions  between  highly  nonlinear 
shock  waves  and  reflecting  boundaries. 

(2)  When  a  high  explosive  of  fixed  mass  at  the  bottom  of  a  pit  is  detonated, 
the  resulting  pressure  wave  has  less  strength  outside  the  pit  than  that 
which  would  have  existed  had  the  blast  not  been  contained  by  the 
pit.  This  was  clearly  shown  in  Fig.  6,  where  one  can  see  that  at  any 
given  distance  from  the  source  of  the  explosion,  the  pressure  at  the 
leading  edge  of  the  blast  wave  for  contained  explosions  is  less  than  the 
leading  edge  pressure  for  an  uncontained  explosion. 

(3)  For  explosions  in  pits  of  variable  radius  and  constant  height,  there  exists 
a  pit  radius  which  will  minimize  the  pressure  at  a  given  distance  from  the 
charge  center.  Figure  6  displayed  the  leading  edge  pressure  versus  pit 
radius  at  two  constant  distances  from  the  charge  center.  It  was  clearly 
demonstrated  that  a  minimum  pressure  does  exist,  which  implies  that 
there  is  a  pit  geometry  for  maximum  effectiveness. 

It  is  clear  that  we  should  look  for  experimental  data  which  could  be  com¬ 
pared  with  our  results  prior  to  attempting  complicated  and  time-consuming 
refinements  to  the  studies  carried  out.  Lacking  such  data,  some  modifica¬ 
tions  could  be  made. 

An  important  improvement  to  our  study  would  be  to  introduce  absorbing 
walls,  rather  than  to  continue  to  use  totally  reflecting  walls.  This  would  have 
several  ramifications.  First,  the  model  would  be  more  realistic.  In  the  actual 
process  of  shock  reflection,  perfect  reflectors  dc  not  exist.  Accurate  model¬ 
ling  of  walls  which  absorb^  a  portion  of  the  shock  wave’s  incident  energy 
could  only  improve  the  calculated  results.  Secondly,  absorbing  walls  would 
most  likely  attenuate  the  blast  wave  in  the  air.  This  is  a  simple  application  of 
the  principle  of  conservation  of  energy.  A  portion  of  the  blast  wave’s  energy 


24 


would  go  into  the  inreveraible  compression  of  the  absorbing  wall.  This  loss  of 
energy  would  weaken  the  blast  wave,  and  thus  increase  pit  effectiveness. 

It  should  be  noted  that  an  appropriate  boundary  condition  for  the  inter¬ 
face  where  the  solid  wall  material  contacts  the  outer  edge  of  the  domain 
would  have  to  be  determined.  No  such  condition  would  be  necessary  at  in¬ 
terior  interfaces  of  wall  material  and  gas;  interior  interfaces  are  handled 
internally  by  CSQ. 

Finally,  it  should  be  noted  that  viscous  effects  in  the  gases  were  assumed 
to  be  negligible.  As  the  wave  progresses  into  the  far-field,  viscous  attenua¬ 
tion  and  vibrational  relaxation  become  the  more  dominant  mechanisms 
for  weakening  the  wave  [8].  It  may  be  profitable  to  examine  these  effects 
at  a  later  date. 
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Summary 

ExperimenU  cited  in  this  peper  reveal  that  -aqueous  foams  are  good  attenuators  of  blast  waves 
and  the  resulting  noise.  A  model  is  preaentuu  which  describes  the  behavior  of  an  explosively  pro¬ 
duced  blest  wave  propagating  through  aqueous  foam.  The  equation  of  state  for  an  air/water  mix¬ 
ture  is  developed  with  specific  attention  to  details  of  liquid  water  compressibility.  Solutions  of  the 
conservation  equations  in  a  spherically  one-dimensional  form  were  performed  using  a  finite-dif¬ 
ference  wave  propagation  code.  Results  are  presented  that  indicate  the  effect  of  the  foam  expan¬ 
sion  ratio  at  well  at  the  dimensionless  foam  depth  on  the  blast  attenuation.  The  (limited) 
comparison  of  decibel  level  attenuation  between  the  model  and  the  experiments  shows  good 
agreement. 


Introduotion 

Spherically  symmetric  blast  waves  resulting  from  explosions  in  air  can  cause 
serious  damage  to  structures  located  many  charge  radii  from  the  center.  In 
addition,  the  blasts  also  produce  significant  levels  of  environmental  noise  at 
distances  beyond  the  region  of  structural  damage.  Consequently,  the  areas 
where  blast-producing  activities  ( such  as  demolition  work  and  ordnance  dis¬ 
posal)  can  be  conducted  safely  are  limited. 

When  a  detonation  wave  propagating  through  a  condensed  explosive  reaches 
the  air/explosive  interface,  an  intense  shock  wave  with  pressures  of  the  order 
of  hundreds  of  atmospheres  is  propagated  radially  outward  through  the  air.  It 
has  been  shown  [  1  ]  that  the  strength  of  the  blast  wave  can  be  greatly  atten¬ 
uated  by  surrounding  the  explosive  charge  with  aqueous  foam.  An  aqueous 
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foam  consists  of  a  matrix  of  thin  sheets  of  water  encapsulating  tiny  pockets  of 
air.  By  adding  a  surfactant  to  the  water  to  increase  the  surface  tension,  foam 
can  be  produced  in  a  wide  variety  of  expansion  ratios.  These  foams  are  most 
commonly  produced  by  commercial  fire-fighting  equipment.  The  expi^ion 
ratio  is  defined  as  the  ratio  of  foam  volume  to  liquid  volume 

a  =  vtM  (1. 

Raspet  and  Griffiths  [  1  ]  summarize  past  experimental  work  on  shock  atten¬ 
uation  and  include  extensive  data  on  far  field  attenuation  of  peak  pressure,  flat 
weighted  integrated  sound  exposure  levels  (FSEL)  and  C-weighted  exposure 
levels  ( CSEL )  for  a  variety  of  charge  masses,  foam  depths  and  foam  expansion 
ratios.  Also  presented  are  scaling  laws  which  relate  the  foam  depth,  foam  den¬ 
sity  and  charge  mass  to  the  noise  reduction  levels. 

Figures  la-c  show  reduction  levels  for  various  depths  of  30 ;  1  expansion  ratio 
foam  [1].  In  these  figures  the  foam  depth  is  scaled  to  the  cube  root  of  the 
charge  mass.  The  experiments  were  conducted  using  three  different  charge 
masses,  0.11  kg,  0.67  kg  and  2.27  kg.  The  reductions  are  plotted  in  dBs  (deci¬ 
bels)  with  the  peak  pressure  level  (fofined  as  20  log  (PJ../P„)  where  P„..  is 
the  mATimiim  pnssure  and  Po  is  a  reference  prassure,  Po= 20  ^a.  The  sound 
exposure  level  (SEL)  is  defined  as 

SEL=101og^~  (2) 

where  to  <  reference  time  defined  to  be  one  second.  The  integral  in  eqn.  ( 2 ) 
is  performed  over  the  entire  duration  of  the  wave  form,  both  positive  and  neg¬ 
ative  phase.  The  FSEL  is  calculated  with  no  jflrequency  weighting  and  the  CSEL 
is  calculated  with  a  standard  C-weighting  filter.  In  all  configurations,  two  trials 
were  performed  and  FSEL,  CSEL  and  peak  pressure  were  measured  at  both 
60  m  and  120  m  from  the  charge  center.  Figures  2a-c  show  the  30 ;  1  data  com¬ 
bined  with  data  for  a  much  less  dense  260;  1  foam.  Here,  the  authors  included 
foam  density  pt  in  the  scaled  depth.  In  Figs.  2a-c  the  dimensionless  foam  depth 
is  defined  as 

X  =  Ar[p,/MV>^  (3) 

where  pt  is  the  foam  density  [kg/m’ I,  dr  is  the  geometrically  averaged  foam 
depth  [m] ,  and  M  is  the  mass  of  the  explosive  [kg]  (TNT  equivalent) . 

In  addition  to  peak  pressure  and  sound  exposure  level  reduction,  Ref.  [1] 
also  notes  a  decrease  in  both  the  total  and  positive  phase  durations  of  the  far 
field  recorded  wave  form  when  foam  is  the  wave  propagation  medium.  Plotted 
as  a  function  of  scaled  foam  depth  these  dB  r^uct’ons  were  20%  and  6%, 
respectively.  As  shown  in  Fig.  1,  the  maximum  reductions  in  CSEL,  FSEL  and 
peak  pressure  were  limited  to  about  10  dB. 

An  interesting  result  presented  in  Ref.  [  1  ]  is  that  the  attenuation  is  shown 
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water  components  do  not  achieve  immediate  equilibrium.  This  leuens  the 
degree  of  attenuation. 

Much  work  has  bean  done  to  explain  the  effeotiveneM  of  foam  aa  a  shock 
attenuator  [2-4] .  These  works  include  theories  based  on  multiple  reflections 
from  bubble  durfaces  and  broa^ning  of  the  shock  due  to  bubble  resonances.  It 
should  1m  noted  that  these  mechanisms  only  occur  in  the  acoustic  or  near 
acoustic  range  of  overpressure  and  are  probably  not  applicable  to  the  extremely 
large  overpreasurea  encountered  in  explosively  initiated  shook  waves.  It  was 
pointed  out  by  Raspet  and  Oriffitha  [1]  that  the  minimum  shock  strength 
needed  to  shatter  the  foam  structure  is  on  the  order  of  several  hundred  kPa. 


Fig.  3.  SchtBUtic  of  tiplotlva/fotm  oonflguration. 

Other  experiments  [6,6]  heve  cited  the  viqporization  of  the  water  and  the 
quenching  of  the  afterburn  as  mechanisms  which  reduce  the  delayed  energy 
imparted  to  the  wave.  We  believe  that  the  energy  contribution  by  afterburn  is 
a  small  percentage  of  the  energy  delivered  by  the  initial  detonation  wave  and 
is  probably  not  the  main  mechanism  for  attenuation.  The  effects  of  vaporiza¬ 
tion  will  be  examined  in  more  detail  later  in  this  study. 

The  work  presented  here  is  by  no  means  a  detailed  study  of  the  pore  surface 
structure  within  the  foam.  Instead,  emphasis  is  placed  on  the  hydrodynamic 
wave  interactions  within  the  medium  and  between  material  interfaces.  Based 
on  the  strength  of  the  explosively  driven  shock  wave  through  the  foam,  the 
hydrodynamic  model  appears  to  be  a  valid  approach.  The  following  sections 
will  outline  the  problem  analysis,  numerical  solution  technique  and  some 
results. 

Analysis 

Figure  3  deHnes  the  problem  at  hand.  Located  at  the  r  s  0  origin  is  a  spherical 
charge  of  explosive  with  radius  Oq.  Surrounding  the  explosive  is  a  foam  with 
expansion  ratio  a  and  depth  dr.  Surrounding  the  foam  is  air.  For  this  analysis, 
the  initial  time  is  taken  as  the  instant  the  detonation  wave  reaches  the  explo- 
sive/foam  interface  (rssOo).  The  form  of  the  detonation  wave  at  this  time  is 
assumed  to  be  the  classical  form  of  Taylor  [7]  for  a  spherical  self-similar  det¬ 
onation  wave  profile.  The  magnitude  of  the  transmitted  and  reflected  waves  at 
the  foam  boundary  are  calculated  from  the  Rankine-Hugoniot  jump  conditions 
and  the  equations  of  state  for  the  two  materials  (explosive,  foam) .  In  order  to 
study  the  wave  motion  beyond  this  initial  time,  one  must  solve  the  equations 
of  motion,  equations  of  state,  >jnd  the  appropriate  boundary  and  initial  con¬ 
ditions  for  the  entire  flow  field.  For  a  one-dimensional,  spherically  symmetric 
analysis  the  equations  of  motion  in  Eulerian  form  are  written  as 
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Coruervation  of  mass 


dp  ,  dp  ,  du  .  „pu  . 

(4) 

Conservation  of  momentum 

du  ,  du  1  dP  . 
dt  dr  p  dr 

(6) 

Conservation  of  energy 

dE^  dE  P  f  dp  ^  dp\  ^ 
dt  ^  dr  p*  \dr  dt)  “ 

(6) 

Equation  of  state 

E^E{P,p) 

(7) 

At  rssO  a  reflected  boundary  condition  is  imposed  and  at  r=ii  a  transmitted 
boundary  condition  is  used.  Internal  boundary  conditions  at  the  material 
interfaces  requite  continuity  of  particle  velocity,  u  and  pressure. 

Along  with  appropriate  initial  conditions,  the  three  nonlinear  partial  differ¬ 
ential  equations  and  the  equation  of  state  form  a  set  of  four  equatioim  to  be 
solved  for  the  four  unknown  variables  E,  P,  u,  and  p.  In  order  to  obtain  an 
analytical  solution  to  this  sat  of  equations,  some  very  limiting  assumptions 
must  be  made  [&-ia].Areviewofmuohoft^workcanbefoundinRef.  [14]. 

The  solution  technique  used  in  the  work  performed  here  is  based  upon  the 
finite  difference  solution  to  the  governing  differential  equations  and  consti¬ 
tutive  relations.  A  brief  discussion  of  the  solution  technique  follows. 

Finite  difference  solution 

The  finite  difference  code  used  to  solve  the  blast  attenuation  processes  is  an 
adaptation  of  the  one-dimensional  WONDY  V  code  developed  at  Sandia  Lab¬ 
oratories  [  16  ] .  A  detailed  description  of  the  operation  of  the  code  can  be  found 
in  Ref.  [16].  A  brief  discussion  of  the  computational  scheme,  stability  and 
ideal  form  of  the  state  equation  will  be  presented  here. 

The  Anite  difference  code  is  used  to  solve  the  set  of  one-dimensional  equa¬ 
tions  of  motion  in  spherical  geometry.  The  governing  equations  written  in 
Lagrangian  form  are 

Conservation  of  mass 

m  =  Wo  (8) 
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(  Conservation  of  momentum 


tf^onaervation  of  energy  * 

Equation  of  state  * 

P  =  P(£,P)  (11) 

Here,  m  represents  the  mass,  p  is  the  density,  9  is  a  viscous  stress,  P  is  the 
pressure,  and  E  is  the  specific  internal  energy. 

In  the  finite  difference  approximation  to  the  differential  equations  all  quan¬ 
tities  sampled  are  material  particles  at  discrete  times.  The  differential  equa¬ 
tions  are  written  in  difference  form  by  the  use  of  centered,  second-order  analogs 
over  a  staggered  conmutational  grid.  The  space  variables,  a,  u,  and  x  (accel¬ 
eration,  velocity,  ana  spatial  location)  are  located  at  the  cell  boundaries,  and 
the  thermodynamic  variables,  P,  E  and  p  are  centered  in  each  cell. 

^Since  the  grid  resolution  of  the  code  cannot  be  made  small  enough  to  accu- 
rhtely  resolve  the  shock  waves  thickness,  an  apparent  viscous  stress  g  is  intro¬ 
duced.  This  prevents  the  wave  form  ^m  overtaking  itself  and  increases 
computational  stability  by  spreading  the  discontinuity  across  several  cells. 
Shock  waves  in  the  finite  difference  solution  are  recognised  as  very  steep  but 
finite  gradients  in  the  solution.  The  form  of  g  used  in  this  work  is 

«  =  C.C.|  +  C.,(i|)’  if|>0  (12.) 


(12a) 


q  =  0 


(12b) 


Here,  Cj  and  C2  are  constants  [  16  ]  and  C.  is  the  local  sound  Velocity.  In  addi¬ 
tion  to  providing  numerical  stability,  the  apparent  viscous  stress  also  satisfies 
the  entropy  production  across  the  shock  front  as  dictated  by  the  second  law  of 
thermodynamics. 


Equation  of  state  for  aqueous  foam 


When  a  shock  passe'j  through  a  liquid-gas  mixture,  the  liquid  requires  a 
finite  time  to  equilibrate  velocity  and  temperature  with  the  gas.  In  the  relax¬ 
ation  zone,  differences  in  velocity  and  temperature  between  the  phases  cause 
momentum  and  heat  transfer  which  can  have  important  effects  on  the  result¬ 
ing  two-phase  flow  field.  Often  these  processes  proceed  very  rapidly,  particu- 
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larly  when  one  phase  is  finely  dispersed  in  th«  other.  When  this  happens  it  can 
be  assumed  that  equilibrium  is  reached  at  the  shock  front  and  that  the  two 
component  systems  can  be  cbnsidered  a  homogeneous  peeudofluid  that  obeys 
the  usual  equations  of  single  component  flow. 

As  stated  earlier,  aqueous  foam  consists  of  a  matrix  of  thin  sheets  of  water 
encapsulating  small  pockets  of  air.  In  this  analysis  it  is  assumed  that  when  a 
strong  shock  hits  these  thin  sheets  they  are  shattered  into  mlcrodroplets  by 
the  viscous  stress  across  the  shock  plane.  This  assumption  is  elaborated  in  Ref. 
[  1  ] .  Because  of  the  droplets’  sIm,  they  equilibrate  very  rapidly  with  the  flowing 
gas.  This  assumption  enables  one  to  consider  the  system  as  an  homogeneous 
entity,  and  thus  use  mass  averaged  thermodynamic  properties  to  describe  the 
system.  These  properties  are  weighted  averages  and  are  not  necessarily  the 
same  as  the  properties  of  either  phase.  Pore  surface  structure  is  not  treated  in 
this  work  bMause  of  the  strong  shock  overpressures. 

In  order  to  analyse  shock  propagation  through  foam  using  a  flnite  difference 
method,  one  must  first  develop  an  equation  of  state  for  the  homogeneous  pseu* 
dofluid  using  average  properties  of  the  air  and  water  components.  The  equa¬ 
tion  is  then  put  in  a  form  most  amenable  to  bydrooode  calculation,  namely  eqn . 
(11). 

The  task  here  is  to  find  the  fiuid  pressure  given  the  fluid  density  and  fluid 
internal  energy  using  both  the  ideal  gas  equation  of  state  with  a  nonconstant 
specific  heat  for  air  and  an  equation  of  state  for  water  of  the  form 
[  P,  H  ]  w/  ( T,  u ) .  The  water  equation  of  state  was  supplied  by  Sandia  Labo¬ 
ratories  [  16  ]  and  was  found  to  be  consistent  to  the  fifth  ilgnifloant  digit  when 
compared  with  tabulated  values  flrom  steam  tables  by  Keenan,  Keyes,  Hill  and 
Moore  (16]. 

The  equilibrium  solution  is  found  by  latisi^ing  the  mixture  mass  (volume) 
and  energy  relations 

mutursencrgy  m  (l>x)fiw  (13b) 

mixtunmcui  Vt  m  xv^,+  il-x)u^  (13b) 

where  x  amass  firaction  of  air,  £  a  specific  Internal  energy,  and  o  a  specific 
volume. 

The  subscripts  f,  a  and  w  stand  for  fluid  mixtures,  air  and  water,  respectively. 
In  addition,  we  have  the  ftinctional  relations 

fi.  «/(r)  -js.(r)  (14a) 

E^^f{T,vJ  mEMuJ  (14b) 

u.  a  P7’/P.  (14c) 
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Fl|.  4.  Shook  ovirpmsur*  for  varloui  •xp«n*ion  ratio  foam*.  B'or  all  oaiai.  Jr/a„  w  40. 


“/(7’iU*)  as  P^(T,y^)  (14d) 

Aiiuining  that  the  preeaure  and  temperature  of  both  phaaee  are  equal,  we  can 
rewrite  the  maaa  and  energy  relatione  ai 

fi,  »  *£?,(T)  + (1 -*)£;»( T,0w)  (ISa) 

RT 

^p^(r,vj  ^^  “ * ^ *^**  ) 

The  two  equation!  are  coupled  through  the  fluid  temperature,  T,  and  the 
water  npeciflc  volume,  u^.  This  provides  two  equations  in  two  unknowns,  T 
and  v^,  The  solution  will  yield  the  fluid  pressure  through  the  water  equation 
of  state,  P*  (T,u^), 

Since  the  internal  energy  of  water  is  dependent  on  both  temperature  and 
apecinc  volume  ( due  to  phase  change )  the  temperature  cannot  be  solved  directly 
from  the  energy  equation.  If  this  were  the  case,  pressure  could  be  solved  for  by 
satisfying  the  mass  relation,  varying  Since  this  is  not  possible,  an  iterative 
scheme  must  be  employed  to  solve  the  two  equations  simultaneously. 

Computed  results 

Figure  4  is  a  plot  of  the  peak  shock  presb.ire  as  a  function  of  distance  from 
the  charge  center  r=0.  Three  different  expansion  ratio  foams  (60,200,360) 
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an  ahown  in  the  figun.  Each  case  illuetrated  in  Fig.  4  has  a  foam  depth  of 
Jr/Oo  a  40.  Figure  4  also  includes  a  plot  of  the  locus  of  the  peak  shock  pressure 
in  an  air  medium  (no  foam  aaoo)  for  the  same  charge  [11].  The  general 
trend  is  for  the  lower  expansion  ratio  foams  to  start  at  a  hig^r  shock  pressure 
than  the  air  and  higher  expansion  foams.  The  higher  initial  shock  pressures 
for  the  lower  expansion  ratio  foams  is  a  result  of  the  foam  shock  impedance 
and  the  product  gas/foam  interface  boundary  condition.  Figure  4  also  shows  a 
break  in  the  rate  of  pressure  decrease  at  about  8  charge  radii  for  the  60 : 1  foam; 
6  charge  radii  for  the  200:1  foam;  and  10  charge  radii  for  the  360: 1  foam.  The 
rate  of  pressure  decay  increases  after  this  point.  Evaluation  of  the  thermody* 
namic  properties  shows  that  this  break  occurs  when  the  lead  shook  is  no  longer 
of  sufficient  strength  to  vaporise  the  foam.  This  attenuation  mechanism 
becomes  evident  when  studying  the  pressure,  energy  and  velocity  histories  of 
the  attenuating  flow  fields. 

Figures  6-7  ^ow  pressure,  energy  and  density  profiles  for  a  constant  expan¬ 
sion  ratio  (2(X):  1)  foam  where  the  foam  depth  is  10  charge  radii.  In  all  six 
different  depths  were  studied  6, 10,  16,  20,  SO  and  40  charge  radii,  with  the 
solutions  integrated  to  a  distance  of  90  charge  radii.  Thermodynamic  profiles 
for  the  remaining  cases  are  not  shown  here,  but  can  be  found  in  I^f.  [14]. 
These  calculations  were  performed  in  order  to  investigate  the  blast  wave 
behavior  at  the  foam/air  interface  and  to  observe  the  effects  of  different  foam 
depths  on  the  far  field  wave  form.  The  initial  blast  wave  for  these  calculations 
wasdeterminedfrom  an  equation  of  state  forTNTused  by  Erode  [13]  starting 
at  the  initial  conditions  specified  by  Taylor  [7] . 

The  results  shown  in  Figs.  6-7  show  the  radial  distance  scaled  to  the  initial 
charge  radius.  Pressure  is  scaled  to  ambient  pressure  ( 0. 101  MPa )  and  density, 
energy  and  time  are  left  unsealed.  Also  t  represents  the  elapsed  time  (/m  )  since 
the  detonation  wave  reached  the  explosive/foam  interface.  It  should  be  noted 
that  prior  to  solving  the  flow  equations  with  foam  as  the  propagation  medium, 
a  test  case  was  run  for  a  standard  charge  in  dry  air.  This  particular  test  of  the 
finite  difference  code  was  made  since  the  corresponding  experimental  results 
are  well  documented  [  17  ] .  The  peak  overpressures  for  this  case  were  in  favor¬ 
able  agreement  with  experimental  observations  in  both  the  near-fieid  and  far- 
field. 

Examining  the  results  shown  in  Figs.  6-7,  one  can  see  that  there  is  no  sig¬ 
nificant  difference  in  the  pressure  profiles  other  than  that  the  second  shock  is 
attenuated  more  and  is  farther  into  the  positive  region  as  the  depth  increases. 
This  does  not  explain  the  change  in  the  rate  of  attenuation.  However,  inspec¬ 
tion  of  the  internal  energy  variations  offers  an  explanation  of  the  phenomena. 

The  increase  in  internal  energy  above  the  ambient  is  the  last  influence  that 
the  blast  wave  produces.  The  kinetic  and  thermal  energy  of  the  initial  explosive 
products  that  goes  into  producing  the  blast  wave  will  be  distributed  in  the  final 
state  as  an  increase  in  the  internal  energy  of  the  surrounding  medium.  When 


Fig.  6.  Shock  overprcHure  hiitory  for  200 : 1  exptnilon  ratio  foam,  dr/tta=  10. 
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Fig.  6.  Specific  internal  energy  hlatory  for  200 : 1  expansion  foam,  Jr/a„  ^10. 


1 


I 


i 


the  shock  wave  passes  a  given  location,  it  leaves  the  material  at  a  pressure  and 
energy  determined  by  the  jump  conditions  and  the  equation  of  state.  When  the 
medium  expands  back  to  ambient  pressure,  it  returns  some  of  this  energy  to 
the  wave  propagating  ahead  of  it  and  retains  some  residual  energy  as  a  result 
of  the  irreversible  heating  caused  by  molecular  shearing  across  the  shock  front. 
As  the  strength  of  the  shock  diminishes,  the  residual  energy  also  diminishes. 


36 


X  LOCATION  r/Qo 

Fig.  7.  Dtniity  hUtory  for  200 : 1  •spanaion  ratio  foam,  Jr/a„ « 10. 


In  this  weak  shock  case  the  wave  can  be  thought  of  as  propagating  acoustically 
with  the  attenuation  mostly  due  to  spherical  divergence.  The  spike  in  the  inter¬ 
nal  energy-distance  plots  (Fig,  7)  is  due  to  the  density  discontinuity  occurring 
at  the  material  interface. 

In  order  to  make  a  comparison  of  the  computed  foam  depth  attenuation 
effectiveness  with  the  data  reported  by  Ref.  [  1  ] ,  the  data  shown  in  Fig.  8  are 
presented  as  reduction  in  dB  peak  level  in  Fig.  9.  Here,  the  dB  reduction  is 
relative  to  the  overpressure  predicted  for  an  r/Oo  of  zero.  Figure  8  shows  this 
to  be  185  ,8  dB.  Thus,  for  example,  a  normalized  foam  depth  of  10  would  give  a 
dB  reduction  of  186.8 -183.2=:  2.6.  The  predicted  dB  reductions  for  five  dif¬ 
ferent  cases  of  r/Oo  are  shown  in  Fig.  9.  Also  shown  (as  a  solid  line)  for  com¬ 
parison  is  the  measured  far-field  peak  level  dB  reduction  for  an  a =250  foam, 
taken  from  Fig.  3  of  Ref.  [  1  ] . 

A  least-squares  fit  of  the  predicted  Peak  Level  Reduction  (PLR)dB  for  nor¬ 
malized  foam  depths  less  than  26  is, 

(PLP.)aB  =  0,28  -0.04  (16) 

In  order  to  compare  this  with  the  data  presented  in  Ref.  [  1  ]  (as  shown  in 
Fig.  9)  it  was  necessary  to  first  carry  out  some  conversions,  since  the  data 
(Fig.  3  of  Ref.  [1])  are  presented  as  (PLR)dB  versus  scaled  foam  depth 
( ) .  This  depth  is  the  ratio  of  foam  depth  divided  by  the  cube  root  of  the 
explosive  charge  mass.  In  terms  of  that  parameter,  the  best-fit  equation  for  the 
a  =  250  foam  was  a  far-field  dfi  reduction  that  is  given  by 
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Fig.  8.  Predicted  "far-field"  overpreesure  decibel  level  as  a  function  of  normalized  foam  thickness 
(a  a  200  foam). 

(PLR)dB  -  6.33  (<dr,/M'/3 -0.03)  forO<dr,/M‘/3  <1.6  (17) 

Since  the  explosive  used  in  the  experiments  was  C-4,  not  TNT,  it  was  neces¬ 
sary  to  scale  the  mass  of  C-4  to  an  equivalent  TNT  mass  of  TNT,  AfTNT=  1-35 
Mc.^.  Therefore,  to  determine  the  equivalent  charge  radius,  ( as  used  in  the 
computer  prediction ) 

=  (Mtnt/1.34)'''''  =  l,73(4;rAhNT/3)‘'=*a„.rNT  (18) 

Thus  1  kg  of  C-4  explosive  is  equivalent  to  a  mass-of  TNT  with  aoTNT=0.04  m. 

When  substituting  this  relation  into  eqn.  (17 )  bne  obtains  the  experimental 
best  least-squares  fit,  as 

(PLR)dB  =  0.366  (Jr, /oo)- 0.19  for  0 r/a„  < 26  (19) 

This  equation  is  shown  by  the  solid  line  in  Fig.  9  and  compares  favorably  with 
eqn.  (16). 

Conclusions 

The  results  presented  show  an  increase  in  the  rate  of  attenuation  after  the 
shock  no  longer  vaporizes  the  water  component  of  the  mixture,  leading  one  to 
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Fig.  9.  Comparison  of  ths  pradictad  dacibsi  paak  lavtl  raduction  (in  tha  “flur-flald”)  with  aspari- 
manta  from  Raf.  ( 1  ]  (solid  lina ) . 

concluda  that  vaporization  can  in  fact  be  detrimental  to  maximum  attenua¬ 
tion.  When  the  liquid  component  of  the  eyatem  vaporizea,  the  air  is  no  longer 
loaded  with  a  relatively  incompressible  material  and  the  shock  speed  increases. 
This  enables  seconded  shocks  to  catch  up  to  and  reinforce  the  main  shock. 
This  aspect  of  phase  change  is  also  detrimental  to  attenuation  in  the  inter¬ 
mediate  field.  The  greater  propagation  speed  of  the  secondary  shocks  may  also 
account  for  the  reduced  duration  of  the  negative  phase. 

The  many  wave  reflections  off  the  foam/air  interface  produce  a  complicated 
waveform  in  this  region.  However,  it  was  noted  that  these  disturbances  rapidly 
decay  into  the  air  region  and  are  small  compared  to  the  peak^disturbance.  For 
smaller  foam  depths  the  transmitted  pressure  is  still  high,  and  nonacoustic 
attenuation  by  the  air  adjacent  to  the  foam  occurs. 

A  significant  pressure  drop  occurs  at  the  foam/air  interface  for  small  foam 
depths,  and  in  this  respect,  the  impedence  mismatch  between  the  air  and  the 
foam  is  important.  However,  for  large  foam  depths,  where  the  attenuation  due 
to  shock  dissipation  occurs  solely  in  this  region,  the  effects  of  reflections  from 
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the  foam/air  surface  and  the  impedence  mismatch  have  little  influence  on  the 
far  field  waveform.  ■ 

Two  factors  are  respo^ble  for  determining  the  amount  of  residual  heat  left 
in  the  medium  once  thlftshock  has  passed.  One  is  the  thermodynamic  proper¬ 
ties  that  influence  the^Bhditions  across  the  shock  discontinuity.  For  maxi- 
mum  attenuation  (for  a. given  pressure)  one  would  wish  to  maximize  the  rise 
in  internal  energy  across  the  shock  and  minimize  the  drop  in  energy  during 
isentropic  expansion- back  to  ambient  pressure.  This  effect  is  determined  solely 
by  the  equation  of  State  of -the  medium'.  The  other  contributing  factor  is  that, 
for  the  same  material,  a'  higher  pressure  shock  produces  a  larger  amount  of 
residual  energy.  ThuC  fQr  large  foam  depths  where  interface  interactions  are 
minimal,  the  fwOinp^  important  factors  contributing  to  high  attenuation  are 
a  maximum  overpressure  and  the  ability  to  produce  a  high  residual 

energy.  In  ed  $;tionj  the  lower  expansion  ratio  foams  attenuate  more  than  light 
foams  for  twq  possible  reasons.  The  first  is  that  the  initial  overpressure  is 
higher  and  the  second"  ia’ that  for  a  given  pressure,  more  residual  energy  is 
produced  It  is’  likely  that  both  factors  contribute. 

Beyonda  certain  foam  depth  the  amount  of  residual  energy  left  in  the  region 
of  expanded  foam  dtbps  sharply.  Most  of  the  energy  imparted  to  the  foam  by 
the  shock  is  then  returned  to  the  v/ave  upon  expansion.  Attenuation  past  this 
point  is  due  mostly  to  divergence. 
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SUMMARY 

This  work  investigates  the  effectiveness  of  several  different  materials  on  the 
reduction  of  blast  noise  produced  by  detonating  high  explosives,  In  each  case, 
the  explosive  charge  is  centered  in  a  cube  of  blast  reducing  material,  The 
materials  .selected  have  quite  different  physical  properties  ( e.g.  steel  wool, 
straw  and  plastic  bubble  pack )  and  hence,  the  results  provide  some  might  into 
the  dominant  mechanism  for  energy  transfer  between  the  shocked  air  and  the 
blast  reducing  materials.  In  all  ca.se.s  reported,  the  flat-weighted  .sound 
expo.sure  level  ( FSEL).  C-weighted .sound expo.sure  level  ( CSEL),  and  the 
peak  .sound  level  ( Pp)  were  measured  at  four  po.sitions,  2  at  38  m  and  2  at 
76  m,  Based  on  the  experimental  data  presented,  all  four  materials  tested 
.scaled  as  a  function  of  the  geometrically  averaged  material  depth,  the 
material  density  and  the  explosive  ma.s.ses  (kg-TNT  equivalent}.  Further¬ 
more.  in  all  ca.ses  studied  the  peak  level  .scales  more  e.xactly  than  do  the  energy 
integrals  CSEL  and  FSEL. 


INTRODUCTION 

It  has  been  recognized  for  some  time  that  the  intense  sound  levels  produced 
by  detonating  high  explosives  can  be  hazardous  for  the  surrounding 

*  To  whom  all  correspondence  should  be  addressed. 
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environment.  With  many  of  the  ordnance  test  and  dispy.<8al  sites  being  slowly 
encroached  upon  by  ever-expanding  nearby  communities,  the  task  of  blast 
noise  reduction  takes  on  growing  importance.  For  the  past  several  years  the 
Acoustics  Team  at  the  U.S.  iWy  Construction  Engineering  Research 
Laboratory  (USACERL)  has  been  investigating  blast  wave  attenuation. 
Much  of  the  work  emphasized  reducing  the  blast  energy  while  it  is  still 
within  several  charge  radii  of  the  source.  One  approach  is  to  surround  the 
explosive  charge  with  an  energy  absorbing  heterogeneous  material  (e.g. 
aqueous  foam,  fibrous  materials,  granular  materials).  In  brief,  a  multiphase 
medium  aids  blast  attenuation  through  increased  energy  dissipation. 

Several  papers  over  the  recent  years*  ~  ^  have  reported  on  the  use  of  both 
high  and  low  expansion  ratio  aqueous  foams  as  blast-reducing 

agents.  The  foam  referred  to  in  these  papers  is  generated  by  means  of  a  batch 
foamer  and  is  similar  in  composition  to  those  used  in  Are  Aghting.  Raspet 
and  GriAAths*  reported  experiments  showing  the  aqueous  foam  acts  as  a 
fairly  good  blast-reducing  agent.  Numerical  calculations  presented  by 
Panezak,  Butler  and  Krier^  illustrate  the  effect  of  the  water  droplet 
vaporization  on  the  overall  attenuation  process.  It  was  concluded  in  Ref.  2 
that  the  post-shock  vaporization  of  the  water  is  actually  detrimental  to 
maximum  attenuation.  This  observation  is  contrary  to  the  speculations 
made  in  previous  work  on  attenuation  of  blast  waves  in  foam.  Evans, 
Jankowski  and  Hirl6man^  suggested  that  for  certain  expansion  ratio  foams, 
mass  loading  is  the  principal  factor  in  attenuation.  More  recent  work  by 
Powers  and  Krier*  reported  on  the  attenuation  effectiveness  of  detonating 
the  explosive  charges  in  an  open  pit  as  opposed  to  the  open  atmosphere.  This 
work  was  carried  out  in  order  to  optimize  the  pit  geometry  for  blast  noise 
reduction. 

The  work  presented  in  this  paper  will  focus  on  the  effect  several  energy¬ 
absorbing  heterogeneous  materials  have  on  the  reduction  of  intermediate 
Aeld  (i.e.  >  ISOdB)  blast  noise.  Test  data  for  several  different  materials  are 
presented.  The  materials  selected  have  quite  different  physical  properties 
(e.g.  steel  wool,  Aberglass,  straw  and  plastic  bubble  pack)  and  hence,  the 
experimental  results  should  provide  some  insight  into  the  dominant 
mechanism  for  energy  transfer  between  the  shocked  air  and  the  blast- 
reducing  materials.  In  all  cases  reported,  the  Aat-weighted  sound  exposure 
level  (FSEL),  C-weighted  sound  exposure  level  [CSEL),  and  the  peak  sound 
level  (PEAK)  were  measured  at  four  positions  in  the  farAeld,  2  at  38  m  and  2 
at  76  m.  The  exact  locations  relative  to  the  explosive  charge  are  shown  in  Fig. 
1.  A  total  of  24  different  conAgurations  were  tested  with  3-4  repetitions  for 
each  conAguration.  The  parameters  varied  in  the  test  were  the  mass  of 
explosive  driver,  the  mass  and  type  of  attenuation  material,  and  the  depth  of 
attenuation  material. 
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Fig.  I.  Microphone  luyout  Tor  test  range  3J,  Fort  Leonard  Wood.  The  microphone 
locutions  arc  indicated  as  IM|,  M],  Mj  and  M4. 


Although  the  data  are  somewhat  limited  in  the  range  of  charge  mass, 
experimental  results  indicate  that  the  noise  reduction  produced  by  the  four 
materials  scales  reasonably  well  with  the  dimensionless  material  depth  X* 
developed  tor  the  study  of  aqueous  foam.*  This  scaling  parameter  {X*) 
includes  attenuation  effects  due  to  the  material  thickness,  apparent  material 
density,  and  mass  of  the  explosive  driver.  In  some  cases  presented,  there  is 
evidence  that  noise  reduction  is  also  dependent  on  fiber  size.  However,  when 
data  from  all  four  materials  are  combined,  the  fiber-size  dependency  appears 
to  be  secondary  when  compared  to  the  mass  loading  effect.  These  data  are 
presented  in  the  following  sections  along  with  a  brief  description  of  the  test 
procedure. 


TEST  PROCEDURE 

All  measurements  were  performed  using  the  same  expermental  layout  (Fig. 
I ).  Endevco  piezo-resistive  microphones  v/ere  used  with  CERL-constructed 
line  drivers;  The  microphones  were  mounted  on  tripods  1-2  m  above  ground 
level.  They  are  indicated  in  Fig.  1  by  M1-M4.  Lines  were  run  from  the 
microphones  to  the  equipment  van  for  recording  and  analysis.  The  C- 
weighted  sound  exposure  level,  the  flat-weighted  sound  exposure  level,  and 
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Pig.  1  Schematic  of  structure  used  to  confine  the  blast  attenuation  materials.  In  all  cases, 

(3ft  or4A). 

the  peak  level  were  measured  with  a  CERL  True  Integrating  Environmental 
Noise  Monitor  and  Sound  Exposure  Level  Meter  for  each  station,  and  the 
signals  recorded  on  an  Ampex  2230 14*track  FM  recorder.  By  definition,  the 
sound  exposure  level  is  given  as 

SEL»101og[f/>*d//Pgro]  (1) 

where  Pq  «ind  to  are  reference  pressure  (20/iPa)  and  time  (1  s),  respectively. 
Peak  pressure  is  defined  in  terms  of  a  decibel  level  (dB) 

PE^^:™201og,o[(Pp“Po)/^]  (2) 

where  is  the  measured  peak  pressure  (Pa).  The  above-mentioned 
instruments  were  calibrated  before  each  test  using  a  Bruel  and  KJaer 
pistonphone  and  the  calibration  was  rechecked  after  the  test. 

A  wood  frame  structure  with  covering  made  of  lightweight  polyethylene 
plastic  sheeting  (Fig.  2)  was  used  to  hold  the  blast-reducing  material  in  place 
prior  to  detonating  the  explosive.  This  arrangement  is  similar  to  the  one  used 
by  Raspet  and  Griffiths'  in  their  previous  work,  in  which  both  high  and  low 
expansion  ratio  aqueous  foams  were  used  as  the  blast-reducing  material.  In 
the  experiments  reported  here,  the  polyethylene  cubes  were  either  0-76S  m^ 
(27  ft’)  or  1-812  m’  (64  ft’)  in  volume.  The  explosive  material  was  C-4,  a 
commonly  used  plastic  explosive  with  an  effective  energy  yield  approxi¬ 
mately  1-36  times  that  of  TNT.  Calculated  detonation  properties  for  C-4  are 
presented  in  Table  1.  In  all  firings,  either  O-S  stick  or  1  stick  (0-57  kg)  of  C-4 

TABLE  I 

Detonalion  Properties  of  C-4  Explosive’ 


Detonation  velocity  (CJ)  8  37  mm  j«s ' ' 

Detonation  pressure  (CJ)  2S-7  CPa 

Bulk  density  of  material  I  66  g  cm  ~  ’ 

Heat  of  detonation  (theoretical)  3-86  MJ  kg' ' 
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was  used  as  the  blast  source.  In  order  to  minimize  ground  eflects,  the 
explosive  charge  was  mounted  on  a  vertical  post  in  the  center  of  the  cubic 
structure.  The  enclosing  structure  was  completely  flilled  with  one  of  the  blast- 
reducing  materials  and  the  charge  was  detonated  remotely.  The  polyethy¬ 
lene  and  wood  enclosure  was  chosen  since  it  is  rigid  enough  to  support  the 
heavier  materials  (ix.  40  kg  of  steel  wool),  yet  does  not  offer  enough 
resistance  to  act  as  an  additional  attenuation  mechanism. 


TEST  RESULTS 

Four  different  blast-reducing  materials  were  chosen  for  this  series  of  tests. 
They  consi.sted  of  steel  wool,  fiberglass,  straw  and  a  plastic  bubble  pack 
commonly  used  for  packaging  and  shipping  delicate  items.  Two  different 
grades  of  the  steel  wool  were  used,  f|i0000  grade  with  a  fiber  diameter  (d)  of 
approximately  ISO/tm,  and  $3  grade  with  an  effective  liber  diameter  of 
240 /um.  The  ffOOOO  grade  has  a  circular  cross-section,  while  the  1^3  grade  is 
more  rectangular.  Both  were  produced  in  roll  form  as  opposed  to  pads.  The 
two  different  grades  will  be  referred  to  as  fine  and  coarse  steel  wool, 
respectively,  in  the  tables  and  figures  to  follow.  Steel  wool  was  chosen  as  a 
potential  blast-reducing  material  for  four  primary  reasons,  the  high  density, 
high  heat  capacity,  high  thermal  conductivity  and  high  surface-to-volume 
ratio  of  the  fibers.  Theoretically,  all  four  of  these  features  should  aid  in  the 
blast  attenuation  through  increased  two-phase  energy  transfer  between  the 
shocked  gas  and  the  solid  material. 

The  fiberglass  was  also  tested  in  two  different  sizes,  a  fine  grade  (d  5iim) 
and  a  coarse  grade  (</»  10/im).  It  has  a  density  about  1/3  of  steel  and  a 
specific  heat  about  twice  that  of  steel,  and  hence  should  have  different 
attenuation  properties.  In  addition,  the  thermal  conductivity  is  significantly 
different  from  that  of  steel  wool.  The  two  different  fiber  diameters  were 
chosen  in  order  to  study  the  effect  of  specific  surface  area  on  blast  reduction. 
In  addition  to  the  steel  wool  and  fiberglass,  one  grade  of  straw  with  an 
external  diameter  of  approximately  2  5  mm  was  tested,  and  the  bubble  pack 
was  tested  for  two  different  bubble  sizes  (</■  14'5mm  and  (/■•290mm). 
Straw  was  selected  since  at  least  one  U.S.  Army  installation  has  used  straw  as 
a  noise  reducing  material.^  It  should  be  noted  that  the  straw  was  cellular  in 
structure  with  an  approximate  internal  diameter  of  50 /im.  Because  of  the 
cellular  structure  of  straw,  it  was  difficult  to  determine  the  exact  surface  area 
exposed  to  the  flow  of  hot  product  ga.ses  from  the  detonation.  In  summary, 
Table  2  lists  the  materials  investigated  along  with  some  of  the  physical 
properties  wc  feel  play  a  key  role  in  analyzing  blast  wave  attenuation. 

Table  3  summarizes  the  test  results.  Listed  for  each  experiment  is  the  muss 
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TABLE  2 

BI«ii>r«<luoin|  Material  Propertiei*'* 


Mimiiil 
ikml/iiUtn  • 

Miiiifrliil 

mmitn'f 

Dluituflrr 

(/Mill 

OtCl.t//,!' 
i*/r»i  *) 

C'mimili'iiy 
(ifm  '/t  'l 

Sfifi'l/ir  Imt 

SihI  wool  mat) 

II 

190 

7750 

«0‘6 

49V 

SiHl  vrool  luoitrM) 

13 

340 

7790 

M)'9 

499 

BubbI*  puck  Umull) 

1.^ 

14900* 

9)0 

0-)H 

14)0 

Bubbit  puck  |lur|«l 

14 

39000* 

9.10 

0'3I 

14.90 

Sifuw 

19 

90 

4) 

00)9 

1  )‘‘0 

Hibtriluu  lllnu) 

lA 

9 

3900 

00)6 

Vbl 

Plbcrvluiu  (cuurMi 

17 

10 

191)0 

00)6 

961 

*  Thu  iilm«n>lun  NprcMinu  iht  ulnglt  hubhli  diiimtitr  Th*  muurUI  hint  u  ihlk'kni)iii  nrO’l  mm. 


TABLE  3 

Teit  Dalu  for  Blail  Wave  Attenuation 


Test 

iwmber 

Material* 

number 

LPEAK 

UtB) 

iiCSEL 

idB) 

^FSEL 

UIB) 

Charged 

IS  sllek) 

Material 
muss  ikg) 

a,  h,  ('t 
(,/» 

20 

SI 

I2'l 

I2'6 

m 

0'5 

47'67 

3-0 

38 

ll^ 
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0'5 
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49 

S7 
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8'9 
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10 
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42 

ss 

I0'9 
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05 

33  14 
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47 

#7 
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07 

«l 

96 
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34 

#6 
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80 
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40 

S2 

90 

93 
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05 
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30 

14 
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86 

82 
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10 
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32 

SI 
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57 
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3042 

4-0 

01 

SI 
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90 

97 

10 

29'51 

30 

to 

S5 

82 

8'2 

90 

05 
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36 

S6 

81 

7'3 

78 

10 

36'77 

4-0 

09 

S2 
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7'9 

99 

05 

29-51 

3-0 

43 

S7 

71 

7'l 

7  1 

0'5 

14-53 

3-0 

12 

S5 

70 

6'4 

6-3 

10 

29-06 

3-0 

04 

S2 

5'9 

63 

70 

10 
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16 

S3 
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4'6 

49 

0'5 

16-80 
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22 

S2 

51 

3'4 

62 
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3-0 

25 

S4 

46 

35 

3'2 

10 

15-44 
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26 

S3 

40 

30 

31 

10 
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S3 
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20 
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28 

S4 
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22 

03 
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- . 

..  -  - 

— 

- — - 

— 

-  - 

•  Sec  Tiible  2  for  mutcriul  type, 
t  One  slick  C-4  -  0'567  5  kg, 
t  Dimension  of  cubic  enclosure. 
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of  explosive  muteriul  used,  muss  of  uttenuution  materiul,  the  material 
designation  (From  Table  2),  and  the  decibel  reductions  recorded  for  peak 
level  {APEAKl  AFSEL  and  ACSEL  The  reductions  were  determined  from 
the  difference  between  the  control  charge  and  the  test  charge  so  the 
variations  due  to  charge  muss,  compositions,  and  temperature  would  be 
minimized.  Where  more  than  one  trial  occurred,  the  maximum  difference  in 
reduction  is  listed.  From  a  safety  standpoint,  the  minimum  value  is  probably 
of  more  iitterest.  The  maximum  variation  usually  occurred  for  the  peak  level 
since  this  measure  is  more  sensitive  to  random  processes  than  the  integrated 
measures  {FSEL  and  CSEL),  On  the  average,  the  reduction  measured  at  the 
closer  microphone  (38  m)  were  one  to  two  dB  greater  than  the  data  measured 
at  the  distant  microphones  (76  m).  This  difference  is  due  to  the  more  rapid 
non-linear  decay  of  the  unmuffled  event.  At  first  glance,  Table  3  illustrates 
several  key  mechanisms  for  material  attenuation  effectiveness.  Material 
muss,  depth  ufmuteriai  and  in  some  cases  muteriul  fiber  diameter  all  appear 
to  govern  the  amount  of  attenuation  a  given  material  produces.  These 
features  become  more  evident  when  the  data  are  categorized  according  to 
muteriul  type.  With  regard  to  a  comparison  between  different  types  of 
materials,  it  appears  that  the  plastic  bubble  pack  is  less  effective  than  the 
more  dense  steel  wool,  straw  and  fiberglass.  These  somewhat  general 
observations  need  to  be  investigated  in  more  detail. 

Figures  3(aHc)  present  the  peak  level,  FSEL  and  CSEL  reductions  for  the 
tests  where  bubble  pack  was  used  as  the  attenuating  material.  The  data  for 
both  small  (f|i3)  and  large  (fji4)  bubble  pack  are  combined  on  these  plots,  and 
although  limited  in  the  number  of  data  points,  they  show  good  agreement 
scaled  as  a  function  of  the  dimensionless  material  depth: 

(3) 

Here,  /  is  the  geometrically  averaged  material  depth: 

l  =  0-5(ahc)''^  (4) 

where  a,  h  and  c  arc  the  enclosure  dimensions  (.see  Fig.  2).  is  the  material 
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bulk  density  (Rgm"-*)  and  is  the  explosive  mass  (kg-TNT  equivalent). 
The  dimensionless  material  depth  given  in  eqn  (3)  was  used  in  previous 
work'  *^  to  scale  the  attenuation  properties  of  both  high  and  low  expansion 
ratio  aqueous  foams.  It  appears  to  be  a  fairly  representative  scaling  law  that 
incorporates  the  material  mass,  depth  and  charge  weight,  all  important 
factors  in  blast  wave  attenuation.  For  bubble  pack,  there  does  not  appear  to 
be  any  significant  effect  of  bubble  size  (i(3  vc  ^4)  on  the  level  of  attenuation. 
Table  4  presents  the  correlation  coefficients  and  linear  regression  formulas 
for  the  three  sets  of  data  shown  in  Figs  3(a)-(c).  The  bubble  pack  shows  a 
peak  level  reduction  of  about  5  dB  per  dimensionless  material  depth,  a  FSEL 
reduction  of  4dB/A'*  and  a  CSEL  reduction  of  4-7  dB/ A'*.  The  previous 
work  with  aqueous  foam*  showed  reductions  of  about  4dB/X*  for  peak, 
FSEL  and  CSEL 


49 


R.  Rtixpci,  R  B.  Butlvr,  F.  Juluini 


filled  with  a  30: 1  expansion  ratio  foam  produced  with  a  Mearl  OT 10  batch 
foamer.  By  definition,'  the  expansion  ratio  is 

0i=K/K  (3) 

where  K,  represents  the  total  volume  occupied  by  the  foam,  and  1^,  is  the 
liquid  volume.  For  future  comparisons,  the  test  configuration  used  by 
Raspet  and  Griffiths'  is  shown  in  Fig.  2.  Here  the  aqueous  foam  is  contained 
in  a  less-confining  plastic  cube  surrounding  the  explosive. 

In  the  present  investigation,  tests  were  performed  with  i  stick  of  explosives 
(0-285  kg)  in  the  0-9 1-m  culvert  with  two  foam  depths  (0-86  and  0-55  m),  and 
with  i  stick  and  1  stick  in  the  1'22-m  diameter  culvert  with  three  foam  depths 
(M6, 0-86  and  0-55  m).  One  stick  of  C-4  has  a  mass  equal  to  0-57  kg  (0*78  kg 
TNT  equivalent).  Table  1  contains  the  results  of  these  tests.  In  addition  to  the 
C-4  mass  and  foam  dimensions  (</-  diameter,  h  =  height),  CSEL  reduction 
{ACSEL),  FSEL  reduction  (AF5£L)  and  peak  noise  level  reduction 
(APEAK)  are  listed  for  the  three  different  combinations  of  charge  size  and 
culvert  size. 

The  CSEL,  FSEL  and  PEAK  reductions  were  determined  from  the 
difference  between  a  control  charge  and  the  test  charge,  so  that  variations 
due  to  charge  mass,  composition  and  temperature  would  be  minimized.  This 
is  the  same  test  procedure  as  used  in  Ref.  1. 

Three  qualitative  features  of  the  blast  wave  reductions  can  be  seen  from 
the  results  in  Table  1.  First,  increased  foam  depth  for  a  given  culvert  and 
explosive  configuration  results  in  increased  noise  reduction.  This  feature  is 
evident  in  the  comparison  of  test  lOia  with  101b,  as  well  as  the  comparison 
of  the  set  including  tests  102a,  102b,  102c,  103a,  103b  and  103c.  A  second 
observation  is  that  an  increased  culvert  diameter  for  a  fixed  depth  results  in 
greater  noise  reduction.  Finally,  increased  charge  mass  results  in  decreased 
noise  reduction.  These  dependences  are  the  same  as  those  found  in  Raspet 
and  Griffiths'  for  unconfined  explosives.  To  display  those  similarities  and 


Fig.  I  Schematic  of  test  conAguration  used  by  Raspet  and  Griffiths.'  Here  the  Foam  is 
confined  by  plastic  sheeting  rather  than  the  rigid  metal  culvert. 
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Reducing  Intermediate  Mast  noise 


Examination  of  the  experimental  data  for  materials  #3  and  #4  suggests 
several  possible  mechanisms  for  the  attenuation  of  the  blast  wave  as  it 
propagates  through  the  plastic  bubble  pack.  The  first  is  based  on  a  thermal 
analysis  of  the  plastic  material.  Since  the  bubble  pack  is  a  combustible 
material  with  a  relatively  high  specific  surface  area,  part  of  the  shock  wave 
energy  can  go  into  heating  up  and  melting  the  plastic.  Subsequent  «t»crgy 
release  from  the  plastic  may  occur  after  the  wave  has  propagated  well 
beyond  the  enclosure  volume.  Another  possibility  is  that  the  bubble  pack 
acts  to  diffract  the  shock  wave  as  it  propagates  through  the  material.  This  is  a 
result  of  the  shock  wave  encountering  the  density  discontinuity  (impedance 
mismatch)  present  in  the  bubble  pack/air  matrix.  The  air  within  each  bubble 
is  at  a  slightly  higher  pressure  than  the  surrounding  air  which  fills  the  voids 
and  hence,  there  exists  a  density  discontinuity.  To  determine  if  the  proposed 
attenuation  mechanisms  are  of  primary  or  secondary  significance,  the 
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attenuation  effectiveness  of  other  materials  should  be  examined.  For  clarity, 
the  data  for  steel  wool  (taken  from  Table  3)  are  given  in  Table  5. 

Figures  4(aHc)  display  the  reductions  for  the  two  different  grades  of  steel 
wool,  fine  ($  1 )  and  coarse  ($2).  It  is  interesting  to  note  that  the  fine  and  coarse 
steel  wool  show  slightly  different  reductions  for  the  same  dimensionless 
material  depth  (A”*).  Taking  into  account  the  definition  of  X*,  this  implies 
that  for  the  same  depth  of  material,  mass  of  material  and  explosive  mass,  the 
two  grades  of  steel  wool  behaved  differently  in  reducing  fhe  blast  noise.  The 
only  difference  between  the  two  is  the  fiber  size,  or  more  specifically  the 
surface-to-volume  ratio  of  the  fibers.  Recall  that  both  ihe  small  and  large 
bubble  pack  scaled  to  X*  and  hence  did  not  show  a  dependency  of 
attenuation  on  a  characteristic  dimension.  This  observation  at  first  leads  one 
to  consider  incorporating  a  surface  area  term  and  possibly  a  material  heat 
capacity  term  in  the  proposed  scaling  law. 
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TABLE  5 

T«st  Data  far  Steel  Wool  (from  Table  3) 
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Fig.  4.  Noise  reductions  for  two  grades  of  steel  wool,  O.  Fine;  y.  coarse.  Scaled  to  A'*,  the 
fine  grade  shows  a  greater  reduction  for  all  metrics,  |u). 
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Fig.  4 — conid  (b) 


With  regard  to  the  thermal  effects,  a  lumped  parameter  heat  transfer 
analysis  yields  a  post*shock  thermal  equilibration  time  of  35  ms  for  the  fine 
fibers  and  86ms  for  the  coa'-se  fibers.*  This  analysis  is  based  on  a  10 MPa 
shock  passing  through  an  air/steel  wool  two-phase  medium.  When 
neighboring  fibers  are  considered,  the  equilibration  times  are  reduced  by  a 
factor  of  about  20.  The  quoted  equilibration  times  are  only  a  first-order 
approximation,  but  illustrate  an  interesting  point.  Because  the  fine  fibers 
equilibrate  with  the  shocked  gases  much  more  rapidly  than  the  coarse  fibers, 
there  is  a  greater  decrease  in  available  energy  to  drive  the  shock  wave 
through  the  remainder  of  the  material.  However,  for  a  strong  shock 
propagating  through  steel  wool,  experimental  data^  shows  that  the  shock 
decay  time  is  of  the  order  of  1  ms.  Since  the  predicted  equilibration  times  are 
of  the  order  of  30-80  ms,  this  implies  that  heat  transfer  effects  are  not 
important  in  dissipating  the  energy.  In  fact,  an  adiabatic  assumption  would 
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suffice.  Note  that  this  analysis  does  not  exclude  surface  area  from  the 
attenuation  model.  Surface  area  viscous  effects  are  still  a  possible 
attenuation  mechanism. 

All  test  data  are  combined  in  Figs  5{aHc).  Here,  the  peak  level,  FSEL  and 
CSEL  reductions  are  again  plotted  as  a  function  of  the  dimensionless 
material  depth  X*.  Based  on  the  data  presented  here,  the  X*  scaling  appears 
to  be  a  reasonable  method  for  combining  data  from  several  materials,  quite 
different  in  physical  properties  (see  Table  2),  loading  densities  and  being 
driven  by  different  charge  weights.  However,  it  is  obvious  that  the  scaling  is 
not  exact  and  other  material  properties  should  be  considered.  Of  the  three 
metrics,  the  peak  level  reduction  appears  to  scale  best  by  the  dimensionless 
X*.  Based  on  the  comb  ned  data  shown  in  Figs  5{aHc),  there  is  a  strong 
indication  that  the  material  mass  loadmg  has  the  strongest  influence  on  blast 
attenuation. 
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Fig.  5.  Noise  reductions  lor  ull  materials  sealed  to  ihc  dimensionless  material  depth.  X*.  (a). 

This,  however,  is  limited  to  the  very  small  range  of  charge  weights 
considered  in  these  experiments.  Further  experiments  lor  greater  charge 
masses  are  needed  in  order  to  extrapolate  these  data  beyond  a  charge  weight 
of  0'57  kg  of  C-4  explosive. 


CONCLUSIONS 

The  work  presented  in  this  paper  represents  an  effort  to  identify  the 
dominant  mechanism  for  blast  wave  attenuation  in  various  heterogeneous 
materials.  In  all  cases  the  energy  source  was  from  detonating  either  0-5  or  1 
stick  of  C-4  explosive.  All  measurements  of  peak  sound  level.  FSEL  and 
CStL  were  made  in  the  intermediate  field  environment  at  either  38  m  or 
76  m  Irom  the  energy  source.  Based  on  the  limited  experimental  data,  the 
lollowing  conclusions  are  drawn. 


Redudng  inlermediaie  blast  noise 
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Fig.  5 — conid  (b), 

First,  all  four  of  the  materials  (steel  wool,  fiberglass,  straw  and  bubble 
pack)  investigated  appeared  to  scale  fairly  well  with  the  dimensionless 
material  depth  X*  defined  in  eqn  (6).  This  finding  is  consistent  with 
previous  data  presented*  for  blast  wave  attenuation  in  aqueous  foams.  It 
should  be  noted,  however,  that  this  scaling  law  should  probably  not  be 
extrapolated  to  higher-charge  masses  without  further  data  to  verify  the 
proposed  law. 

The  parameter  X*  is  a  measure  of  the  material  mass  per  charge  mass  and 
appears  to  be  the  best  single  scaling  parameter  available  for  these  types  of 
heterogeneous  materials.  This  is  not  to  imply  that  other  factors  such  as  fiber 
size  and  thermal  properties  of  the  material  do  not  influence  attenuation. 
Materials  that  are  homogeneous  or  have  very  low  void  fractions  would 
probably  not  scale  according  to  this  law.  In  all  cases  studied,  the  peak  level 
scales  more  exactly  than  do  the  energy  integrals  CSEL  and  FSFL,  Although 
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Fig.  5 — contd.  (c). 
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the  characteristic  fiber  size  (d)  influenced  the  degree  of  attenuation  in  some 
cases,  it  was  not  included  in  the  scaling  law.  The  difference  in  specific  surface 
urea  between  the  fiberglass  and  the  steel  wool  is  several  orders  of  magnitude 
and.  as  u  single  parameter,  does  not  describe  the  relative  attenuation 
difference  between  samples.  If  the  fiber  characteristic  length  is  included  in 
future  modeling  efforts,  it  should  be  in  the  form  of  a  dimensionless  viscous 
drag  term.  In  this  study  it  was  observed  that  the  effect  of  fiber  size  is  more 
pronounced  for  tiie  steel  wools  than  it  is  for  the  much  finer  fiberglass.  One 
possible  explanation  is  that  the  material  fiber  size  mav  reach  a  saturation 
level  where  viscous  effects  arc  essentially  independent  of  characteristic 
length. 

In  conclusion,  the  data  presented  provide  the  basis  for  engineering 
estimates  of  blast  wave  noise  reductions  at  distances  between  38  and  76  m 
from  the  source  and  for  charge  weights  up  to  0  57  kg. 
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SUMMARY 

Experiments  were  performed  in  order  to  quantify  the  additional  attenuation 
provided  hy  enclosing  a  blast  reducing  material  ( aqueous  foam )  in  a  rigid 
vessel  ( cylindrical  metal  culvert),  open  at  one  end  to  the  atmosphere.  The 
results  are  compared  with  previoti.dy  reported  data  on  aqueous  foam  where 
the  culvert  was  not  used.  A  total  of  eight  configurations  were  investigated. 
Tests  performed  with  a  0-9 1-m  culvert  section  and  j  stick  (0-285  kg)  ofC-l 
explosive  and  with  a  l-22-m  culvert  section  with  f  stick  and  I  stick  ofC-4.  A 
third  parameter  varied  in  the  trials  was  the  amount  of  foam  used  ( depth  in 
culvert ).  Data  are  presented  for  FSEL,  and  CSEL  and  peak  level  reduction 
.scaled  according  to  a  modified  scaled  foam  depth  dependent  on  the  charge 
weight,  and  height  and  depth  of  the  culvert.  This  modified  scaling  law 
illustrates  the  relative  effectiveness  of  enclosure  depth  and  width  on' the  noise 
reduction.  t  ^ 


1  INTRODUCTION 

The  research  discussed  in  this  paper  deals  primarily  with  the  reduction  of 
blast  wave  overpressures  resulting  from  detonating  high  explosives.  As 
reported  in  previous  work,*'*  potentially  dangerous  sound  levels  can  be 

.Applied  Acoustics  0(X)3-682X/87  S03-50  ■<•  Elsevier  Applied  Science  Publishers  Ltd. 
England.  1987.  Printed  in  Great  Britain  Used  with  permission. 
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mitigated  through  the  use  of  energy  absorbing  materials  at  the  blast  source. 
As  a  continuation  of  the  earlier  efforts  to  understand  the  reduction  of 
intermediate  (ix.  >lS0dB)  noise  levels  produced  by  detonating  high 
explosives,  the  US  Army  Construction  Engineering  Research  Laboratory 
(USACERL)  has  investigated  the  attenuation  levels  produced  by  surround¬ 
ing  an  explosive  charge  with  aqueous  foam  which  is  confined  within  a  rigid 
cylindrical  vessel.  This  research  is  similar  to  the  work  reported  in  Ref.  1  in 
which  the  sound  absorbing  material  (aqueous  foam)  was  supported  by  thin 
plastic  sheeting,  an  enclosure  design  which  presumably  added  no  additional 
attenuation. 

The  experiments  reported  herein  were  performed  in  order  to  quantify  the 
additional  attenuation  provided  by  enclosing  the  blast  reducing  material  (in 
this  case  aqueous  foam)  in  a  rigid  vessel,  open  at  only  one  end  to  the 
atmosphere.  The  results  are  compared  with  the  previously  reported  data  on 
aqueous  foam.‘  Scaling  laws  wilt  illustrate  the  relative  effectiveness  of 
enclosure  depth  and  width  on  the  level  of  noise  reduction.  As  in  the  previous 
worki^'^  the  range  of  charge  weights  is  limited  to  0'28kg<C,„<0‘S7kg. 

2  EXPERIMENTAL  SET-UP 

The  tests  discussed  in  this  paper  were  conducted  at  the  Fort  Leonard  Wood, 
Missouri,  demolitions  training  range.  The  physical  layout  of  the  test  facility 
is  similar  to  the  one  described  by  Raspet,  Butler  and  Jahani.^  As  in  the 
previous  research,  Endevco  piezo-resistive  microphones  were  mounted  on 
tripods  1-2  m  above  ground  level  at  various  distances  from  the  blast  source. 
In  all  four  microphones  were  used,  two  at  38  m  from  the  charge  and  two  at 
76  m.  Each  pair  was  separated  by  90°  relative  to  the  blast  source.  In  each  test 
case,  three  metrics  were  measured  and  recorded  by  the  remote  data 
acquisition  system.  The  C-weighted  sound  exposure  level  {CSEL),  the  flat- 
weighted  sound  exposure  level  (FSEL)  and  the  peak  level  (PEAK)  were 
measured  for  each  of  the  four  stations,  and  the  signals  recorded  on  an 
Ampex  2230 14-track  FM  recorder.  To  assure  reliable  results,  the  system  was 
calibrated  prior  to  and  after  each  test  using  a  Bruel  &  Kjaer  piston-phone. 
By  definition,  the  peak  sound  pressure  level  (PEAK)  and  sound  exposure 
level  (SEL)  are  given  as 

PEAK(6B)  =  20log,o[(/>p-  Po)//’o]  (D 

SELidB)  *=  10  log, 0 Jp"  dt/Plt^  (2) 

where  P^  is  the  peak  thermodynamic  pressure  (Pa),  P^  is  a  reference  pressure 
(Po  =  20/<Pa),  and  Iq  is  a  reference  time  ((„  =  1  s). 
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CHARQE 


Pig.  1.  Schematic  of  test  configuration  used  in  this  work.  The  explosive  charge  is  placed  in  a 
metal  culvert  filled  with  30:1  expansion  ratio  aqueous  foam. 


All  test  results  described  in  the  folibwing  section  (see  Table  1)  were 
performed  with  aqueous  foam  supported  by  a  metal  culvert  section.  This 
particular  configuration  is  shown  in  Fig.  1.  Two  charges  were  used  in  each 
test;  a  control  charge  set  on  a  0‘9-m  high  crushable  post,  and  a  test  charge  set 
at  the  center  of  the  cylinder  on  a  4-cm  high  post.  The  culvert  sections  were 
sunk  into  the  ground  about  6  cm  to  reduce  the  propagation  of  noise  under 
the  cylinder.  To  some  extent,  this  phenomenon  did  occur  since  the  culvert 
section  was  driven  into  the  air  by  the  blast.  A  total  of  eight  different 
configurations  were  investigated.  This  included  tests  that  were  performed 
with  the  0-9  Um  culvert  section  and  ^  stick  (0*285  kg)  of  C-4,  and  with  a 
1-22-m  culvert  section  with  ^  stick  {0*285  kg)  and  1  stick  (0*57  kg)  of  C-4.  The 
third  parameter  varied  in  the  trials  was  the  amount  of  foam.  In  all  cases,  the 
experimenters  attempted  to  keep  the  foam  expansion  ratio,  a,  fixed  at  a  value 
of  30:1  (total  voIumeMiqUid  volume). 

The  material  referred  to  as  C-4  is  a  military  explosive  (91%  RDX  +  6% 
TNT  -I-  3%  other)  with  a  detonation  pressure  Po  =»  25*7  GPa,  detonation 
velocity  0^  =  8*4  mm  ns '  ‘  and  theoretical  heat  of  detonation  of  = 
5*86  MJ  kg“ For  comparison,  it  is  about  1*36  times  as  energetic  as  TNT, 


3  TEST  RESULTS 

As  discussed  earlier,  the  series  of  tests  were  performed  in  order  to  study  the 
effect  of  external  confinement  on  the  blast  reduction  produced  by  aqueous 
foams.  A  schematic  of  the  test  configuration  is  shown  in  Fig.  1 .  The  confining 
vessel  is  a  cylindrical  metal  culvert,  open  at  the  top  end.  During  all  tests,  the 
explosive  charge  was  first  set  in  the  test  cylinder  and  then  the  cylinder  was 
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TABLE  I 

Test  Data  and  Scaled  Foam  Depths 


Test 

number 

Explosive 

mass 

(kg) 

(/ 

(m) 

h 

(mf 

isPEAK 

aCSEL 

(tlB) 

EFSEL 

(dB) 

X 

lOla 

0'28 

0'92 

0'86 

20-6 

18'S 

t7'8 

0*71 

3'75 

lOlb 

a28 

0'92 

0'55 

I6'8 

t4'4 

140 

0-61 

2'77 

102u 

0'28 

l'22 

1  16 

24'9 

21 '6 

19'4 

0*95 

4.99 

102b 

028 

l'22 

0'86 

23'3 

22'2 

20'8 

0‘86 

4'09 

102c 

0'28 

I '22 

0'55 

180 

t5'8 

IS'I 

0'74 

3-04 

103a  . 

057 

t'22 

116 

22-0 

20'4 

19'8 

0-7,5 

3-96 

i03b 

0'57 

1'22 

0'86 

21-1 

19'8 

19-2 

0-68 

3-23 

103c 

0'57 

1-22 

O'SS 

161 

14-4 

14'3 

059 

2-41 

any  diffetences  in  a  quantitative  manner,  it  is  necessary  to  apply  scaling  laws 
to  these  data. 


4  DISCUSSION  OF  RESULTS 

These  tests  were  performed  so  that  the  scaling  laws  developed  in  Raspet  and 
Griffiths'  could  be  applied  to  cases  where  the  foam  is  supported  by  rigid 
walls.  Data  presented  in  Ref.  1,  where  unconfined  aqueous  foam  was  u.sed  as 
the  blast  reducing  agent,  indicated  that  a  scaled  foam  depth 

(4a) 

scaled  foams  of  the  same  density,  and  that  a  dimensionless  foam  depth 
defined  as 

(4b) 

scaled  foams  of  different  densities  quite  well.  In  cqns  (4),  p,  represents  the 
foam  density,  I  represents  the  geometrically  averaged  pit  depth  (see  Fig.  2), 
and  C„  is  the  mass  of  explosive  in  equivalent  kilograms  of  TNT.  This  is  also 
the  same  form  of  the  scaling  used  in  Ref.  3,  where  fiberglass,  steel  wool  and 
other  materials  were  used  in  place  of  foam.  Since  all  the  experiments 
reported  in  Refs  1  and  3  were  for  the  cubic  arrangement  shown  in  Fig.  2,  the 
geometrically  averaged  foam  depth  was  given  as 

/ «  ifl  (3) 

where  a  is  the  dimension  of  the  cube  shown  in  Fig.  2. 

Because  the  current  experimental  layout  is  somewhat  similar  to  that 
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reported  in  Ref's  1  and  3,  and  because  the  observed  trends  are  also  similar,  a 
scaling  law  such  as  eqn  (4b)  is  a  prime  candidate  for  these  new  data. 

Before  this  can  be  accomplished,  a  reasonable  way  of  adjusting  for 
charges  at  the  bottom  of  the  foam  volume,  rather  than  centered  in  the 
volume,  must  be  developed.  Theoretically,  there  are  two  counteracting 
effects.  First,  the  energy  from  the  charge  propagates  into  a  solid  angle  of  In 
rather  than  4n,  thus  increasing  the  effective  charge  by  a  factor  of  two. 
Second,  energy  is  absorbed  by  the  ground,  which  results  in  a  reduction  of  the 
effective  charge  weight.  We  can  use  the  results  of  two  experimental  tests  to 
guide  us  in  choosing  an  effective  weight.  Tests  at  Fort  Leonard  Wood  in 
cor^unction  with  earlier  tests ^  indicated  that  charges  bn  the  ground  were 
quieter  by  3dB  FSEL 

Measurements  of  the  sound  level  reduction  from  charges  set  in  the  culvert 
with  no  blast  noise  reducing  materials  are  given  in  Table  2.  It  is  interesting  to 


TABLE  2 

Reductions  with  No  Foum  in  Culvert 


Metric 

d 

(tn) 

Reductitm 

(dB) 

PEAK 

()‘9| 

-03 

ESEL 

0-91 

to 

CSEL 

0-91 

to 

PEAK 

l'22 

1-3 

FSEL 

1'22 

2'5 

CSEL 

l'22 

2'5 

note  that  the  system  of  culvert  plus  ground  has  a  larger  reduction  on  the 
energy  measures  than  on  the  peak  levels.  This  indicates  that  the  dissipation 
in  this  case  is  possibly  due  to  multiple  reflections  after  the  initial  shock  front 
develops.  In  view  of  the  variation  in  data  and  the  sparsity  of  data,  the  best 
adjustment  is  the  simplest.  We  will  use  the  geometrically  averaged  foam 
depth  divided  by  the  cube  root  of  the  charge  weight.  This  is  equivalent  to 
assuming  that  half  the  energy  of  the  charge  Is  propagated  into  the  ground. 
We  will  scale  the  data  using  this  simple  model  and  note  any  variations  which 
might  be  caused  by  this  particular  choice. 

For  comparison  with  the  previous  data,  the  scaled  foam  depths  {^)  for 
each  of  the  current  tests  are  listed  in  Table  1.  The  CSEL  reduction  versus 
scaled  loam  depth  is  displayed  in  Fig.  3u,  the  FSEL  reduction  in  Fig.  3b  and 
the  peak  reduction  in  Fig.  3c.  Also  displayed  on  these  figures  are  the  lines 
fitted  to  the  uncondned  data  in  Ref.  1. 

For  the  cylindrical  geometry  of  the  culvert,  the  characteristic  foum  depth  / 
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is  given  us  the  equivalent  radius  of  a  sphere  with  the  same  volume  as  the 
cylinder.  That  is. 

'•*  (6) 

Obviously,  these  tests  are  not  extensive  enough  to  provide  detailed  scaling 
laws  giving  reduction  as  a  function  of  diameter,  depth  and  charge  size. 
However,  these  tests  do  provide  an  engineering  estimate  of  the  reduction 
produced  by  confined  charges. 

A  notable  feature  of  this  data  is  the  tendency  of  the  dB  reduction  to 
saturate  for  the  larger  depths.  There  is  a  large  change  in  reduction  from  0  55 
to  0-86  m  depth,  and  a  smaller  change  in  reduction  from  0*86  to  I- 16  m  depth. 
Since  the  foam  is  uncontined  in  the  vertical  directions,  perhaps  vertical 
saturation  is  beginning  to  occur  at  these  depths. 

It  appears  that  equivalent  depth  variations  produce  larger  changes  than 
equivalent  diameter  variations.  That  is,  if  the  foam  volume  is  kept  constant 
so  that  the  scaled  foam  depth  is  constant  but  the  depth  increased  and 
diameter  decreased,  the  reduction  will  increase.  A  more  accurate  theory  may 
need  to  divide  the  volume  dependence  of  scaling  into  an  area  dependence 
and  a  depth  dependence. 

Figures  4a,  b  and  c  show  the  FSEL  reduction,  CSEL  reduction  and  peak 
level  reduction  plotted  us  a  function  of  a  modified  scaled  foam  depth 
defined  us 


where  A  is  the  surface  urea  of  foam  and  h  is  the  depth  of  foam.  The  scaled 
data  shows  good  agreement,  with  the  exception  of  one  erroneous  data  point 
(103b).  When  one  excludes  this  data  point  along  with  the  data  point  that 
appears  to  be  in  the  saturation  region,  a  linear  regression  analysis  gives 
reasonable  results.  These  arc  displayed  in  Table  3  for  all  three  metrics.  Here 
the  reductions  are  417,  405  and  4'54dB/sculed  distance  for  PEAK, 
FSEL  and  CSEL,  respectively. 


table;  .1 

Linear  Rcgrcsilun  Culvulutiuns  fur  Experi- 
mcntul  Dula  Loss  I02u  and  lO.tb:  V-  AX ■¥  B 


Variahk  K 

Correlation 

I'oefficlent 

A 

B 

^PEAK 

(>■99 

417 

5' 5.1 

A  FSEL 

096 

405 

3'44 

\CSEL 

096 

454 

2'4.1 
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Fig.  4. 


s 


<c) 

(c)  Pcuk  level  reductions  for  aqueous  foam  in  a  metal  culvert. 


5  CONCLUSIONS 

In  this  work,  we  have  investigated  the  attenuation  eiTectiveness  of  aqueous 
foam  confined  in  a  rigid  metal  structure.  When  compared  to  previous  data* 
for  unconfined  foam,  it  appears  that  the  rigid  confinement  contributes  to  the 
attenuation  process.  It  may  be  that  the  wave  reflections  off  the  internal  walls, 
as  well  as  the  heterogeneous  medium,  dissipate  energy.  Previous  data' 
showed  that,  for  the  unconfined  case,  the  dimensionless  material  (foam) 
depth  provided  a  reasonabie  scaling  law.  For  the  confined  situation,  it  was 
shown  that  for  several  different  configurations  the  best  scaling  was  with 
foam  depth  raised  to  the  j  power  and  surface  area  to  the  ^  power  in  addition 
to  the  regular  charge  mass  to  the  --i  power.  The  system  of  culvert  plus 
ground  has  a  larger  reduction  in  the  energy  levels  than  in  the  peak  levels. 
This  indicates  that  the  dissipation  in  this  case  is  due  to  multiple  reflections 
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after  the  initial  shock  front  develops.  This  is  consistent  with  the  wave 
reflection  statement  made  above. 

It  should  also  be  noted  that  the  geometrical  arrangement  provided  in 
hese  tests  results  in  a  focusing  of  the  blast  energy  in  the  vertical  direction. 
The  degree  of  attenuation  provided  at  distances  greater  than  76  m  can  only 
be  determined  by  measurements,  since  it  is  not  clear  if  the  readings  at  76  m 
are  affected  by  this  energy. 
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